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ABSTRACT 


Using  a  statistical  regression  model,  a  relationship  was 


found  between  the  area  of  solar  active  regions  (bright  in 

\  2 

hydrogen-alpha,  X *  6562.8  A  in  the  visible)  and  their  10.7 
cm  radio  emission.  Data  used  in  this  analysis  were  the 
brightness-area  histograms  obtained  from  the  Solar  Observing 
Optical  Network  (SOON)  and  daily  drift  scans  obtained  from  the 
Algonquin  Radio  Observatory,  Ottawa,  Canada.  This  data,  for 
the  period  1  May  through  16  July  1978,  was  used  to  develop  the 
following  relationship: 


Algonquin  Scan  Peak  Value  * 

.88782838  +  .00025237604  (Calibrated  SOON  Area) 


with  a  correlation  coefficient,  r  ,  ofx .88509*  Future  models 
for  predicting  the  10.7  cm  solar  flux  are  theorized  based  on 
the  established  relationship. 
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CORRELATING  RADIO  EMISSIONS  VITH  VISIBLE  EMISSIONS  FROM 


INDIVIDUAL  SOLAR  ACTIVE  REGIONS 


I.  Introduction 


Background 

A  definite  correlation  exists  between  hydrogen-alpha  ac¬ 
tive  regions  on  the  sun  and  the  10.7  cm  solar  radio  flux  ob¬ 
served  at  1700Z  daily  at  the  Algonquin  Radio  Observatory,  Ot¬ 
tawa,  Canada  (Refs  7:454  and  8:126).  Since  this  correlation 
has  been  shown,  officers  at  the  Solar  Forecast  Center,  Air 
Force  Global  Weather  Central  (AFGVC),  Offut  Air  Force  Base, 
Nebraska,  want  to  predict  the  actual  10.7  cm  flux  from  obser¬ 
vations  'of  active  regions  with  the  Solar  Observing  Optical 

\ 

Network  (SOON).  Indeed,  one  of  the  ancillary  purposes  of  this 
research  will  be  to  determine  the  feasibility  of  using  SOON 
data  as  an  input  to  a  predictive  model.  The  primary  purpose 
of  this  research,  however,  will  be  to  find  a  relationship 
between  the  Algonquin  drift  scans  of  the  sun  and  the  data  from 
the  SOON  network.  The  SOON  network  has  only  existed  since 
1977,  and  to  date  little  of  the  data  has  been  used  for  any 
research  purposes  (Ref  38).  The  underlying  rational  for  this 
study  is  found  in  an  article  by  Covington,  Legg,  and  Bell 
describing  the  relationship  between  hydrogen-alpha  filtergrams 
-■*  e»  west  drift  scans  of  the  sun  (Ref  9:471-473)*  Instead 


of  using  hydrogen-alpha  filtergrams,  this  data  is  replaced 
vith  the  SOON  data  which  is  also  centered  on  the  hydrogen- 
alpha  wavelength  (  x  M  6562.8  A  ). 


SOON  vs  Algonquin  Scan  Daily  SOON  Visible  Model 

Correlation  Measurements  (Sunspots) 

(This  Study)  I 


Future  10.7  CM  Solar  Flux  Value 


10.7  CM  Solar  Flux,  UV, 
and  X-ray  Correlation 


Future  UV  and  X-ray  Values 

I 

Environmental  Forecasts 
(Atmospheric /Radar /Radio) 


Figure  1.  Prediction  Hierarchy 


Knowing  active  region  area  apd  sunspot  activity, 
forecasters  want  to  develop  an  accurate  way  of  predicting  the 
10.7  cm  flux  from  these  features  when  they  reappear  after 
rotating  out  of  the  field  of  view  of  the  earth.  This  involves 
roughly  a  13-day  lead-time  depending  on  the  region’s  latitude. 
Such  a  time  development  model  of  visible  active  regions  would 
rely  heavily  on  the  proposed  relationship  being  investigated 
by  this  study  as  well  as  current  models  of  sunspot  movement, 
development,  and  decay.  These  components  could  be  synthesized 
into  a  predictive  routine  for  the  10.7  cm  solar  flux  with 
roughly  a  13-day  lead-time. 
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Prediction  of  the  10.7  cm  flux  ia  not  the  end  product; 
rather,  it  is  used  as  an  indicator  of  solar  ultraviolet  (UV) 
and  X-ray  radiation.  There  is  a  documented  correlation 
between  the  10.7  cm  radio  flux,  solar  UV  and  solar  X-ray 
radiation  (Ref  11:2).  The  10.7  cm  flux  is  used  as  a  predicter 
for  these  latter  two  radiations  because  10.7  cm  data  can  be 
collected  on  the  ground  whereas  UV  and  X-ray  measurements  must 
be  done  above  the  atmosphere.  Being  able  to  predict  the  10.7 
cm  flux  in  advance  would  allow  similar  lead-time  prediction  of 
UV  and  X-ray  wavelengths.  Of  primary  interest  are  the  effects 
of  UV  and  X-ray  radiation  on  the  earth's  upper  atmosphere. 
These  include  the  density  of  the  upper  atmosphere  which  deter¬ 
mines  the  drag  experienced  by  low  altitude  satellites.  Also, 
UV  and  X-ray  radiation  affects  the  structure  and  conditions  of 
the  ionosphere.  Many  agencies  are  very  concerned  with 
ionospheric  effects  because  of  the  way  they  influence  radar 
and  radio  transmissions.  Figure  1  highlights  these  inter¬ 
relationships. 

Now,  with  these  factors  in  mind,  a  relationship  must  be 
developed  between  the  SOON  data  and  Algonquin  Radio  Obser¬ 
vatory  drift  scans  of  the  sun'  which  can  be  used  in  future 
studies  to  predict  the  10.7  cm  flux  using  SOON  data.  This 
involves  exploring  growth  and  decay  trends  expected  in  the 
size  of  bright  active  regions  from  one  day  to  the  next,  how 
the  rotation  of  the  sun  affectB  apparent  area  of  these  regions 
at  the  earth  due  to  foreshortening,  viewing  position  (accoun¬ 
ting  for  inclination  of  the  earth's  orbital  plane  and  tilt  of 
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the  solar  axis),  v  itudinal  affects  on  the  sun  itself,  and  if 
the  active  region  is  above  or  below  the  solar  equator.  The 
above  considerations  can  be  gleaned  from  brightness-area 
histograms  from  the  SOON  network.  The  histograms  are  recorded 
daily  for  selected  active  regions  on  the  sun.  All  of  the 
above  factors  can  be  roughly  accounted  for  based  on  a  com¬ 
posite  of  the  individual  histograms. 

Statement  of  the  Problem 

Data  concerning  the  relationship  between  hydrogen-alpha 
active  regions  and  the  total  solar  10.7  cm  radio  flux  has  ex¬ 
isted  since  1977,  but  no  work  has  been  done  to  relate  these 

two  data  bases  in  the  short  term  on  a  region  by  region  basis. 
Most  studies,  -none  involving  SOON  data,  deal  with  longer 
periods  such  as  a  month  and  the  11  year  solar  cycle. 

Ob  j ec tive' of  the  Research 

The  objective  of  this  research  will  be  to  develop  a  rela¬ 
tionship  between  SOON  data  and  Algonquin  drift  scans  which  can 

be  used  as  the  basis  for  a  forecast  model  for  the  10.7  cm 
radio  flux  using  brightness-area  histogram  scans  from  the  SOON 
network.  Basic  hypotheses  involved  concern  whether  or  not  the 
SOON  data  is  related  to  the  Algonquin  scans  and  whether  it  is 
significant  enough  to  be  used  as  a  predictive  tool. 

Methodology 

The  objective  of  this  study  will  be  accomplished  by 
analyzing  data  comprised  of  SOON  brightness-area  histograms 
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and  10*7  cm  radio  flux  data  from  the  National  Oceanic  and  At¬ 
mospheric  Administration  (NOAA),  Boulder,  Colorado,  and  Algon¬ 
quin  Radio  Observatory,  respectively.  The  histograms  and 
solar  scans,  though  taken  at  different  locations,  are  done  at 
roughly  the  same  time  and  can  therefore  be  compared  directly. 
A  regression  model  will  be  developed  with  the  data.  The  model 
will  be  run  with  the  help  of  SPSS  (Statistical  Package  for  the 
Social  Sciences)  regression  packages  and  the  CYBER  1700  com¬ 
puter. 

Limi ta  tions 

As  with  any  new  system,  SOON  has  some  inherent  limita¬ 
tions.  In  the  preliminary  stages  of  this  research,  it  was 
hoped  that  digitized  pictures  of  the  solar  disk  would  be 
available  for  analysis.  This,  however,  was  not  the  case.  The 
sun  is  divided  into  regions  for  SOON  observation,  and  all  are 
not  looked  at  daily.  In  fact,  the  SOON  data  available  is 
limited  to  active  regions  on  any  given  day.  This  means  that 
the  entire  solar  surface  will  not  be  available  for  analysis 
every  day,  and  that  a  composite  of  available  histogram  data 
will  have  to  be  used  to  formulate  a  crude  solar  image  for  com¬ 
parison  with  the  solar  scan  data.  Additionally,  there  are 
some  gaps  in  the  data  due  to  cloudy  weather  conditions  at  the 
SOON  and  Algonquin  sites  and  days  when  observations  could  not 
be  made  because  of  technical  difficulties  (Refs  4  and  30). 
For  active  regions  near  the  limbs  of  the  sun,  large  errors  are 
introduced  into  the  histogram  data  if  the  dark  background  sky 
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is  included  in  the  measurements.  Data  is  generally  recorded 
for  regions  as  they  move  across  the  entire  solar  disk  but  is 
most  accurate  for  those  regions  lying  65°E  and  65°W  of  the 
central  meridian  (Ref  4).  Therefore,  most  of  the  data  from 
.the  solar  limbs  will  have  to  be  eliminated. 

Vith  these  limitations  in  mind,  the  research  was  further 
constrained  by  the  availability  of  the  SOON  data  itself. 
Since  the  system  is  relatively  new,  little  use  has  been  made 
of  the  data.  Because  of  this,  most  of  the  data  is  in  raw 
form,  requiring  much  time  and  effort  to  reduce  to  a  usable 
form  (brightness-area  histograms).  To  date,  the  only 
available  SOON  data  in  reduced  form  covers  the  period  from  May 
1978  to  December  1978.  This  data  base  must  be  further  reduced 
since  only  data  corresponding  to  the  1700Z  Algonquin  drift 
scans  is  relevant.  Therefore,  of  the  available  data  taken 

throughout  the  day,  only  the  readings  (histograms)  as  close  to 

\ 

1700Z  as  possible  for  each  region  observed,  will  be  used  in 
this  study.  Unfortunately,  there  is  one  further  problem. 
Mr.  Arnold  Star,  Detachment  4,  3rd  Weather  Wing,  Holloman  AFB, 
one  of  the  observer/collectors  of  the  SOON  data,  informed  me 
that  he  had  serious  doubts  as  to  the  validity  of  the  data  for 
this  period  in  its  present  histogram  form.  Mr.  Star  feels 
that  the  data  for  the  May  1 978-December  1978  period  may  not  be 
standardized.  As  he  explained,  a  videometer  is  used  to  break 
the  image  into  64  shades  of  gray.  The  pixels  in  each  band  are 
then  counted  and  equated  to  an  area.  The  greatest  number  of 
counts  should  be  centered  on  the  background  level  to  obtain  a 
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normalised  quiet  sun  level.  Early  readings  (including  the  Nay 
1 978-December  1978  period)  were  not  all  normalized  to  the  same 
quiet  sun  level.  Some  were  normalized  to  10  and  some  to  12 
(on  a  0-64  scale).  Mr.  Star  suggested  that  the  data  be  fur¬ 
ther  normalized  to  account  for  all  of  these  effects  (Ref  39). 
This  exercise  is  beyond  the  scope  of  this  study.  Instead,  the 

data  will  be  used  as  is.  More  accurate  correlations  and 

% 

predictions  will  have  to  be  the  result  of  future  studies  in 
this  area  using  later  data  that  is  more  standardized. 

One  other  obvious  limitation  with  this  technique  is  the 
time  involved  in  getting  the  data  into  a  usable  format.  For 
short-term  predictions,  the  speed  with  which  the  data  is 
processed  and  made  available  for  analysis  will  have  to  be 
greatly  increased.  Otherwise,  the  SOON  data  will  be  of  little 
practical  value  for  short  term  predictions. 

Thetfe  are  no  readily  apparent  problems  with  the  10.7  cm 

\ 

solar  flux  data  from  the  Algonquin  Radio  Observatory.  This 
data  has  been  used  as  a  basis  for  prediction  of  the  10.7  cm 
solar  flux  since  1947.  It  has  been  the  standard  basis  for  all 
models  with  earth-based  parameters  (Ref  11:2). 
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II .  Historical  Survey 


SOPH  Da  ta 

The  Solar  Observing  Optical  Network  is  one  of  the  newest 
means  by  which  solar  related  data  can  be  collected  on 
hydrogen-alpha  active  regions.  This  worldwide  system  was 
developed  and  is  operated  by  the  United  States  Air  Force 
(USAF)  Air  Veather  Service.  As  of  now,  there  are  four  sites 
designed  to  collect  data.  Data  collection  began  at  Holloman 
AFB,  New  Mexico,  in  November  1977  and  was  soon  supplemented  by 
data  from  Palehua,  Oahu,  Hawaii,  December  1977,  Ramey,  Puerto 
Rico,  February  1979,  and  Learmonth,  Australia  in  1979.  A 
fifth  site  has  been  proposed  for  somewhere  in  the  Middle  East. 


A  OPERATIONAL  V  PROPOSED  □  UNOER  CONSTRUCTION 
*operoted  by  USAP  Air  Neolher  Service 


Figure  2.  Solar  Observing  Optical  Network  (Ref  4) 
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The  observing  schedule  of  active  regions  at  each  observatory 


is  decided  upon  by  the  duty  observer,  with  suggestions  from 


the  parent  site  at  Holloman,  depending  on  the  level  and  type 


of  activity  and  the  needs  of  real-time  services;  this  deter¬ 


mines  the  order  of  data  on  magnetic  tape  and  companion 


photographic  film. 


As  the  system  is  so  new,  extensive  work  is  still  being 


done  to  determine  how  best  to  use  all  of  the  data.  As  uses 


are  developed,  more  and  more  of  the  data  will  be  processed 


into  the  final  outputs,  which  consist  of  brightness-area 


histograms  (see  Figure  3)»  flare  time-histories,  as  well  as 


REGION  N0:1092  .  McMATU:15266  YY  MM  DD  IIH  MM  SS  Ha  OFFSET 
LAT"  2028  LON  =  7.79  78  5  1  19  58  56  0.0000 


BRIGHTNESS -AREA  PLOT  FROM  HOLLOMAN 


Figure  3»  Brightness-area  scan  of  an  active  region  (Ref  30 1 1  ) 


9 


photographs  of  the  corresponding  solar  phenomena.  The 
brightness-area  histogram  depicted  in  Figure  3  is  a  graphical 
representation  of  all  the  data  provided  in  each  SOON  histogram 
record  set.  Each  record  set  could  be  used  to  construct  a 
histogram  in  this  form.  The  system  is  calibrated  such  that  if 
the  quiet  sun  mode  occurs  in  bin  10  (bin  equals  0-64),  bin  20 
is  twice  as  intense  and  bin  30  is  three  times  as  intense. 
Borrowing  from  the  SOON  descriptive  pamphlets 


The  basic  SOON  telescope  is  a  25-cm  evacuated  refractor 
mounted  on  a  polar  axis  using  a  monochromatic  filter 
(0.5A  bandpass)  centered  on  the  Hydrogen-alpha  absorption 
line  of  the  solar  spectrum.  A  built-in  minicomputer 
totally  directs  the  collection  of  all  data. 


The  digital  H-alpha  brightness-area  data  are  the  major 
data  product  from  a  SOON  telescope.  These  data  present 
brightness  versus  area  information  from  particular  flare- 
producing  regions  on  the  sun.  The  brightness  ranges  from 
dark' shades,  indicating  relatively  cool  absorbing  plasma, 
to  very  bright  shades  indica  ting>.  ho  t  emitting  plasma, 
i.e.,  flares.  Sixty-four  shades  of  brightness  are 
graphically  presented  on  the  horizontal  axis  of  the 
brightness-area  plot.  The  height  of  a  brightness  column 
represents  the  uncorrected  area  of  a  given  brightness 
shade  summed  over  the  entire  region.  The  observer  sets 
threshold  values  for  filament,  quiet  sun,  plage,  and 
flare  according  to  observing  conditions,  and  the  cor¬ 
responding  areas  in  millionths  of  the  solar  disk  (not 
corrected  for  foreshortening)  are  calculated  and  recorded 
on  this  plot.  The  automated  capability  of  the  SOON 
telescope  system  allows  brightness-area  information  to  be 
collected  from  as  many  as  six  regions  per  5-sec  interval. 
The  system  is  usually  programed  to  collect  1  to  3 
brightness-area  scans  per  minute  from  each  selected  re¬ 
gion  (Ref  30) . 


The  10.7  CM  Solar  Radio  Flux 

Radio  emissions  from  the  sun  range  from  8  mm  to  1J  ■  in 
wavelength  of  which  the  10.7  cm  wavelength  (2800  MHz)  is  only 
a  part  (Ref  3:189)*  This  band  is  essentially  blackbody  radia¬ 
tion  from  high  temperature  regions  and  varies  slowly  between 
fairly  narrow  limits  (Ref  28:30).  Also,  part  of  the  10.7  cm 
radio  flux  is  emitted  by  spiraling  electrons  in  the  sun's  mag¬ 
netic  field  and  is  sensitive  to  the  strength  of  the  magnetic 
field  so  that  active  regions  with  strong  magnetic  fields  may 
emit  stronger  in  the  10.7  cm  radio  flux  than  older  less  active 
regions  (Ref  13:21).  These  radio  emissions  are  made  up  of 
three  components; 

a  fairly  constant  basic  component  (quiet  sun  emission), 
occasional  bursts  of  different  types  and  durations, 
closely  associated  with  solar  flares,  and  a  slowly 
varying  component. .. exhibi ting  roughly  a  27-day 
recurrence  period  (Ref  10:734). 

The  basic  component  can  be  found  by 

plotting  the  apparent  temperature  vs  projected  sunspot 
area  and  extrapolating  to  zero  sunspot  area... the  slowly 
varying  component  arises  from  regions  at  a  height  of 
100,000  km  above  the  photosphere  and  [is]  associated  with 
chromospheric  faculae  (Refs  3:189  and  35:223)* 

This  has  been  verified  since,  as  the  sun  rotates  and  varied 
size  spots  and  spot  groups  move  across  the  solar  disk,  the 
slowly  varying  component  varies  accordingly.  Thus,  the  slowly 
varying  component  shows  a  close  association  with  sunspot  num¬ 
ber  and  exhibits  the  27  day  variation  (Ref  35:225). 

Basically  then,  research  indicates  that,  over  and  above 
the  background  basic  component,  the  slowly  varying  component 
of  the  10.7  cm  radio  flux  is  emitted  from  the  corona  and  has  a 
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10.5  year  period  associated  with  it  (Refs  1 1 : 1 ,8; 6 : 31 ;8: 1 29 ; 
and  28:33)-  Often,  the  three  components  described  above  are 
not  broken  out  in  this  manner.  Only  two  components  are 
discussed,  active  and  background  (Ref  1:3232). 

With  these  regions  and  components  in  mind,  a  brief 
discussion  will  be  made  of  the  methods  used  and  values  ob¬ 
tained  for  the  solar  flux.  The  antennas  used  to  obtain  the 
radio  flux  are  of  low  resolution  so  that  the  antenna  beam 
covers  the  whole  solar  disk  and  the  values  obtained  represent 
an  average  over  the  solar  disk  (Ref  15:4)-  These  values  have 
ranged  from  67  x  1 0“22Watts/m2/Hz  to  380  x  1 0“22Wa t ts/m2/Hz 
(Vatts/m  /Hz=Watts  m  Hz  =Jansky)  since  1947  (Ref  36:Section 
5»l)  never  falling  below  65  x  1 0”^Watts/m^/Hz  (Refs  11:4  and 

8:127).  Richards  has  noted  that  "the  values  may  differ  by  as 

22  2 

much  as  100  [x  10“  Wa tts/in  /Hz]  from  one  day  to  the  next  near 

22 

their  maximum,  and  generally  by  no  more  than  6  or  10  [x  10“ 
Watts/m2/Hz]  near  their  minimum  values”  (Ref  36:Section  5»l)- 

The  1700Z  radio  measurements  for  this  study  were  obtained 
from  the  Algonquin  Radio  Observatory,  Ottawa,  Canada,  with  the 
cooperation  of  Mr*  Aunri  Gagnon.  The  drift  curves  are  taken 
with  a  multi-element  interferometer  capable  of  providing  mul¬ 
tiple  fan  beams,  1.5*  of  arc  east-west  by  2°  north-south,  suf¬ 
ficiently  separated  so  that  only  one  response  falls  on  the  sun 
at  any  one  time.  When  this  is  used  in  conjunction  with  an 
offset  grating  of  four  elements,  a  compound  interferometer 
results  capable  of  achieving  a  single-lobed  fan  beam  with  15" 
arc  east-west  (Ref  9:465-466).  This  resolution  produces  scans 
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like  the  one  depicted  in  Appendix  B  (Figure  7)  which  is 
representative  of  all  the  Algonquin  data.  The  peaks  represent 
radio  emissive  regions  associated  with  plages  and  active 
regions . 

10.7  CM  Solar  Radio  Flux  Uses 

The  10.7  cm  solar  radio  flux  data,  having  the  advantages 
of  being  measured  on  the  earth's  surface  and  for  the  longest 
period,  are  being  used  in  many  ways.  Because  of  its  long 
history  of  careful  observation  (Ref  19:327)  and  its  use  as  a 
general  predictor  of  solar  radiation  in  Extreme  Ultraviolet 
(EUV)  and  X-ray  wavelengths  (Refs  1 9 : 327 ; 1 8 : 1 0 ;  and  11:2),  the 
10.7  cm  flux  has  evolved  as  a  well-entrenched  parameter  that 
is  used  by  a  large  number  of  radio  communicators  and  upper 
atmospheric  modelers  for  a  variety  of  operational  applications 
(Refs  1 9  : 327 ; 1 4 : 1 -2 ;  and  8:131)* 

Usefulness .  Descriptions  of  the  usefulness  of  this 
parameter  range  from  it  being  an  excellent  predictor  of  solar 
activity  (Ref  28:33)  to  an  opinion  by  one  expert  that  its  time 
has  passed  (Refs  12  and  13:21  )  and  new  and  better  methods 
should  be  adopted  for  more  accurate  prediction.  Richards, 
commenting  on  the  10.7  cm  flux,  want  as  far  as  to  say  that 
"the  yearly  mean  is  all  that  can  be  reasonably  predicted  since 
the  shorter  period  fluctuations  are  quite  irregular  and  follow 
no  pattern"  (Ref  3&:Section  5,0* 

Pred ic ted  Values .  There  are  similar  feelings  concerning 
the  spectrum  that  the  10.7  cm  flux  is  used  to  predict.  In 
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regard  to  EUV,  many  scientists  feel  that  the  ten-centimeter 
flux  is 

too  gross  a  parameter  to  accurately  evaluate  the  state  of 
the  ionosphere  and  upper  atmosphere  and  that  more 
specific  measurements  in  the  EUV  would  improve  the 
operational  models.  However,  the  lack  of  a  consistent 
and  on-going  EUV  observational  program  and  the  expense 
and  disruption  of  redoing  present  operational  models  (Ref 
19:327) 

precludes  eliminating  the  ten-centimeter  flux  as  a  model  in¬ 
put.  Until  the  advent  of  routine  exoa tmospher ic  EUV  and  X-ray 
satellite  measurements,  the  10.7  cm  flux  was  seen  as  a  fairly 
accurate  predictor  of  these  events  (Ref  28:33)*  Now,  however, 
since  such  measurements  are  more  routine,  numerous  studies 
point  out  that  the  10.7  cm  flux  is  not  a  very  accurate  predic¬ 
tor  at  all  (Refs  1 9 : 327 ; 1 : 3231 ; 1 2 ;  and  13:21)  and  should  not 
be  relied  upon  for  accurate  ionospheric,  upper  atmospheric, 
and  satellite  drag  models.  Certain  correlations  are  very  good 
such  as  the  one  between  the  10.7  cm  flux  and  atmospheric  den¬ 
sity  changes  resulting  from  heating  (Ref  1:3231)  but  this  is 
an  aggregate  effect  over  a  relatively  long  period.  Whenever 
10.7  cm  flux  data  is  smoothed,  it  generally  gives  a  fairly 
good  correlation  in  the  long  term  but  short  term  correlations 
are  very  poor  (Refs  23:1;13:22).  In  this  regard,  there  are 
documented  effects  due  to  solar  rotation  and  whether  measure¬ 
ments  are  taken  near  solar  maximum  or  minimum  (Ref  1:3231)* 

When  dealing  specifically  with  UV,  Donnelly  points  out 

that: 

sunspot  number  and  10*7  cm  solar  radio  flux  are  useful  as 
a  general  guide  to  the  level  of  solar  activity.  However, 
they  are  not  precise  indications  of  the  solar  UV  flux. 
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They  can  give  very  erroneous  UV  estimates  when  they  are 
used  for  particular  days  rather  than  for  average  effects 
(Ref  15:21). 

He  goes  on  to  say  that 

the  10.7  cm  radio  flux  and  sunspot  number  are  of  course 
correlated  with  the  solar  UV  flux  because  all  three  are 
generally  related  to  the  number  and  size  of  active 
regions.  The  more  smoothing  applied  to  the  data,  the 
higher  the  correlation  because  the  differences  due  to  CMD 
[Central  Meridian  Distance]  and  magnetic  field  orienta¬ 
tion  and  strength  distribution  among  active  regions 
become  smoothed  out  to  average  values  (Ref  15:22). 

Such  correlations  are  also  found  when  10.7  cm  flux  and  flares 
are  considered.  Mohler  has  shown  that  the  10.7  cm  daily  flux 
provides  a  general  guide  to  flare  activity  but  that  the  ac¬ 
curacy  varies  for  different  flux  levels  (Ref  29:2).  The  same 
general  finding  applies  to  the  correlation  between  the  quiet 
sun  soft  X-ray  flux  and  the  10.7  cm  radio  flux  as  investigated 
by  Gibson  and  Van  Allen  (Ref  18:8-10).  Finally,  Bossy  and 
Nicolet  use  the  10.7  cm  flux  as  an  input  for  their  predictive 
model  of  Lyman-alpha  irradiance  (Ref  2:907). 

From  the  foregoing  survey,  it  can  be  surmised  that  better 
predictors  are  available  with  documented  validity.  They  will 
eventually  work  their  way  into  various  predictive  routines. 
But,  until  their  data  bases  are  significant  and  continuously 
recorded,  as  the  Ottawa  10.7  cm  flux  data,  their  usefulness 
will  be  limited. 

10.7  CM  Solar  Radio  Flux  Predic  tive  Models 

A  result  of  the  relative  newness  and  lack  of  processing 
of  the  SOON  data,  coupled  with  numerous  conversations  with 
officials  directly  involved  with  collecting  and  processing  the 


15 


data,  it  is  clear  that  no  one  has  attempted  to  develop  a 
prediction  scheme  for  the  10.7  cm  solar  flux  using  the  SOON 
data  (Refs  4;17;  and  22).  There  have,  however,  been  numerous 
other  predictive  techniques  developed  to  predict  the  10.7  cm 
flux  in  the  long  (up  to  a  year),  medium  (month),  and  short 
term  (days).  By  using  SOON  data,  forecasters  hoped  to 
eliminate  some  of  the  observer  dependent  subjectivity  as¬ 
sociated  with  existing  10.7  cm  prediction  methods. 

Cautions .  The  relationship  between  the  Zurich  relative 
sunspot  number,  plage  index,  and  the  10.7  cm  flux  has  been 
established  and  used  by  many  scientists  to  predict  the  10.7  cm 
flux  (Ref  7:454  and  8:126).  A  brief  review  of  efforts  from 
the  long  to  short  term  follows  some  appropriate  cautions  from 
Yu  I.  Vitinskii,  an  expert  in  solar  activity  forecasting. 
Vitinskii  points  out  that  any  prediction  based  on  sunspots  and 
sunspot  areas  suffers  from  inherent  inaccuracy  due  to 
visability  conditions  and  limb  effects.  Atmospheric  condi¬ 
tions  and  proximity  of  spots  to  the  solar  limbs  add  a  great 
deal  of  uncertainty  to  the  number  and  area  of  spots.  If  these 
parameters  are  used  in  prediction  models,  the  output  is  less 
than  optimum  since  errors  are  compounded.  This  is  one  reason 
monthly  and  yearly  values  are  preferred  so  the  values  can  be 
smoothed  to  reduce  inaccuracies  (Ref  42:5)*  Vitinskii's  ob¬ 
servation  related  to  these  long  term  techniques  is  often 
quoted:  "we  have  shown  that  the  reliability  of  the  results 
obtained  using  these  methods  still  leaves  much  to  be  desired" 
(Ref  42:103). 
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Long  Term  Techniques .  The  extensive  effort  of  Slutz  et 
al  is  an  example  of  using  smoothed  Zurich  13  month  mean  sun¬ 
spot  numbers  to  predict  the  10.7  cm  solar  flux.  This  study 
was  supposed  to  be  an  aid  in  solar  activity  prediction  by 
first  developing  a  model  to  predict  sunspot  numbers  several 
years  in  advance  and  then  use  these  values  to  predict  the  cor¬ 
responding  10.7  cm  flux  value.  An  integral  part  of  the  study 
was  therefore  to  determine  which  of  several  mathematical  ex¬ 
pressions  best  described  the  relationship  between  the  10.7  cm 
flux  and  the  Zurich  13  month  mean  sunspot  number.  The  results 
indicated  that  a  second  order  equation  was  acceptable  to 
represent  the  relationship  over  the  12  month  period  since,  up 
to  the  time  of  the  study,  little  consideration  had  been  given 
to  this  relationship  for  periods  less  than  one-half  of  a  solar 
cycle  (l  cycle  is  approximately  11  years)  (Ref  37:87).  The 
long-term'  technique  advocated  in  this  study  was  considered 
valid  enough  to  be  used  as  the  10.7  cm  flux  prediction  method 
in  NASA's  Solar  Activity  Pr edic tion  of  Sunspot  Numbers  (Ref 
24:1)  and  its  subsequent  verification  (Ref  32).  Note, 
however,  that  this  model  and  its  applications  are  used  in  the 
long  term  for  the  13  month  smoothed  10.7  cm  flux  value  and  do 
not  impact  on  the  present  study. 

Some  atmospheric  density  and  satellite  drag  modelers  have 

also  evolved  their  own  methods  for  predicting  the  10.7  cm  flux 

as  well  as  other  parameters  needed  in  their  models.  Richards 

developed  a  piece-wise  linear  model  for  different  levels  of 

the  smoothed  sunspot  number  (R  ): 
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where  F^ ^  ^  is  the  yearly  mean  flux  value  (Ref  36:Section 
5,2).  Slutz'  colleagues,  Stewart  and  Leften,  compared  the 
results  from  their  final  second  order  equation  for  predicting 
the  yearly  mean  10.7  cm  flux  value  (C^g) 

*12  -  63.7451  +  0.728015R12  ♦  .000890443(1^  2)2  (4) 

with  Richard's  equations  values  and  found  the  results  to  be 
very  similar  (Ref  37:107)*  Again,  the  long  term  nature  of 
these  studies  negate  their  usefulness  for  short  term  predic¬ 
tions  • 

As  discussed  by  Gibson  and  Van  Allen,  Xffftus  •  scrmined 

two  components  of  the  10.7  cm  flux  in  his  booK  Radio  As  xronomy 

\ 

(Ref  26:Chapter  3)*  He  showed  that  the  radio  flux  density  at 
earth  from  an  optically  thin  emitting  region  is 

P  =  2.7  x  10"48T_1/2S  (5) 

where  P  is  the  radio  flux  density  at  10.7  cm  at  1  Astronomical 
Unit  in  units  of  10~2^Watts  (m2Hz)“^,  T  is  the  electron  tem- 

perature  in  10°°Kf  and  S  =  / dV  ,  the  volumetric  emission 

-3  /  -3 

measure  in  cm  (N  is  the  electron  number  density  in  cm  and 

© 

dV  is  the  element  of  volume  in  cm^)  (Ref  18:15).  The  volume 
integral  is  taken  over  the  emitting  region,  assumed  to  be  at 
the  same  temperature  T  throughout.  For  an  optically  \,nick 
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emitting  region 


P  -  1 . 1  x  10"19AT  (6) 

2 

where  A  is  the  projected  area  in  cm  of  the  outer  'surface  of 
the  emitting  region  as  viewed  from  earth  (Ref  18:16).  These 
equations  were  developed  using  data  from  a  limited  period 
where  P  ranged  from  90. 1<.P<_261  .7  .  It  seems  as  though  these 

equations  could  be  used  in  the  short  term  depending  on 
variable  availability,  but  this  was  never  done  on  a  regular 
basis . 

Medium  Term  Techniques .  Jack  Slowly  at  the  Smithsonian 

As trophysical  Observatory  developed  the  following  relationship 

for  use  by  NOAA  to  predict  the  monthly  flux  based  on  the 

smoothed  Zurich  relative  sunspot  number  13  month  average. 

This  technique  uses  a  least-square  regression  model  to 

generate  the  predicted  smoothed  R^  values  which  are  then  used 

as  inputs  to  linear  regression  equations  to  compute  Fin  / 

iU.Tvi?/ 

solar  flux  values  (Refs  19:332  and  16:380) 

^10.7d3)  *  49,4  +  °*97R13  ♦  17.6exp(-0.035R13)  •  (7) 

This  equation  was  used  by  Holland,  Rhodes,  and  Euler  to 
extend  values  of  the  10.7  cm  flux  back  to  1749  in  order  to 
develop  a  predictive  model  of  solar  activity  (Ref  23:2). 

Short  Term  Techniques .  For  short  term  predictions, 
forecasters  at  the  Global  Weather  Central  initially  used  a 
regression  technique  to  establish  a  trend  based  on  the  last  3 
solar  rotations  (weighted  toward  the  most  recent  rotation)  and 
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the  trend  for  the  past  two  or  three  days  of  the  current 
rotation.  Then,  using  the  east-west  10.7  cm  scans  from  the 
National  Research  Council,  Ottawa,  to  establish  the  relative 
location  of  radio  plages  and  thereby  estimate  the  effects  of 
solar  rotation,  they  varied  their  prediction  based  on 
experience  (Ref  19:341  and  22).  This  procedure  may  lack  rigor 
but  until  that  time  provided  the  only  very  short  term  forecast 
of  10.7  cm  flux  (see  Appendix  A  for  sample  output)  (Ref 
28:48) . 

The  most  rigorous  models  used  for  predicting  smoothed 
values  for  the  10.7  cm  flux  and  geomagnetic  indices  are  SUN¬ 
SPOT  developed  by  the  NASA  Marshall  Space  Flight  Center,  the 
Air  Force  Global  Weather  Central  27-day  regression  forecast 
subroutine,  and  EST10  developed  by  NOAA.  SUNSPOT  is  used  for 
statistical  estimation  of  future  expected  values  of  the 
monthly  1.0.7  cm  solar  flux  and  the  5  month  smoothed  flux 
values  (^q  j)  (Ref  41  )  •  The  AFGWC  27-day  regression  forecast 
subroutine  predicts  the  10.7  cm  solar  flux  for  the  27  days 
following  the  current  day’s  10.7  cm  flux  (Ref  22).  EST10  is 
used  to  estimate  the  10.7  cm  flux  value  for  1,  2,  and  3  days 
into  the  future.  It  must  be  continually  updated  and  run  in 
order  for  the  short  term  forecasts  to  be  accurate.  To  date, 
this  is  the  most  refined  short  term  predictive  tool.  In  fact, 
it  virtually  ignores  any  long  term  effects  and  produces  very 
subjective  results  since  the  operator  can  change  the  inputs 
based  on  his  experience  (Ref  22).  As  with  each  of  the  other 
models,  these  procedures  do  not  rely  on  any  SOON  data  inputs. 
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III .  Methodology 


Data  Reduction 

Before  the  data  could  be  used,  it  had  to  be  reduced  into 
formats  for  comparison.  The  Algonquin  data,  in  graphic  form, 
of  solar  scans,  could  be  used  for  direct  comparison  with  the 
SOON  data.  The  scans  are  taken  of  the  sun  from  the  west  limb 
to  the  east  over  a  period  of  about  two  to  four  minutes  daily 
at  1700Z  (Ref  9:465  and  27:749)*  See  Appendix  B  for  an  exam¬ 
ple.  The  scans  are  of  the  apparent  sun  without  correcting  for 
the  solar  tilt  in  relation  to  the  ecliptic  (P  angle)  or  the 

heliocentric  location  of  the  sun's  center  (B  and  L  ).  These 

o  o 

angles  did  have  to  be  considered  when  trying  to  make  a  direct 

comparison  with  the  SOON  data  since  P,  B  ,  and  L  are  accoun- 

ted  for  in  the  region  identification  of  this  data.  For  the 

purpose  of  this  analysis,  P,  B  ,  and  L  were  considered  con- 

o  o 

stant  for  the  entire  day.  In  the  ephemeris  they  are  for  0000Z 
whereas  they  are  used  for  corrections  to  measurements  taken  at 
1700Z  (Ref  40:350-553)*  Small  variations  produced  by  this 
difference  were  not  considered  significant. 

Each  brightness-area  histogram  from  the  SOON  network  has 
a  specific  format  (See  Appendix  C).  The  first  card  of  a  four 
or  eight  card  sequence  (single  or  double  length  histogram, 
respectively)  has  all  the  identifying  information  for  the 
particular  region  being  examined  and  calibration  data  for 
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being  able  to  relate  each  histogram  record  to  the  others* 

Since,  on  a  given  day,  in  the  time  frame  of  interest  (near 

1700Z),  the  SOON  instrument  can  sample  numerous  regions  many 

times  each,  the  optimum  approach  was  to  average  the  data  for 

each  region  observed*  The  SOON  program  was  written  for  this 

purpose.  This  program  not  only  averages  the  daily  readings 

near  1700Z  for  each  region  observed,  but  also  corrects  the 

% 

region  area  to  account  for  foreshortening  (See  Appendix  D)  and 
performs  a  coordinate  transformation  on  the  original  SOON 
latitude  and  longitude  to  get  the  data  points  in  the  reference 
frame  of  the  Algonquin  data  (Ref  25:8-15)*  See  Appendix  E  for 
a  derivation  of  the  transformation  algorithm  used.  This 
transformation  allows  for  direct  comparison  of  the  SOON  bright 
regions  with  the  bright  region  peaks  on  the  Algonquin  scans 
and  provides  a  rough  idea  of  the  location  of  the  region  in 

terms  of  solar  radii  distance  from  the  central  meridian.  This 

\ 

is  done  with  a  spherical  to  rectangular  coordinate  transforma- 
tion  (Ref  34:35-38). 

In  order  to  validate  the  SOON  program,  a  sample  set  of 
data  was  selected  to  exercise  the  various  aspects  of  the 
program.  As  with  the  real  data,  it  allowed  for  summing  the 
areas  above  the  quiet  sun  level  for  a  rough  approximation  of 
activity  in  each  region  measured,  selecting  the  appropriate 
daily  P,  Bq ,  and  LQ  readings  from  the  ephemeris,  deciding 
whether  a  single  or  double  length  histogram  was  being  read 
(column  1  of  card  10),  reading  the  region  identifier,  day, 
time,  latitude,  longitude,  wavelength  offset,  quiet  sun  level, 
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and  calibration  constant.  Most  of  these  features  are  then 
used  in  the  program  to  obtain  summary  data  or  as  checks  to 
determine  if  the  data  set  should  be  discarded  (for  example,  a 
wavelength  offset  other  than  zero  causes  the  histogram  to  be 
purged;  wavelength  offset  gives  an  indication  of  how  far  from 
the  hydrogen-alpha  centerline  the  readings  are.  If  it  is  not 
zero,  additional  corrections  must  be  made  to  the  readings 
which  were  beyond  the  scope  of  this  study).  After  accounting 
for  all  of  these  factors,  the  data  was  in  a  form  for  direct 
comparison . 

Da ta  Comparison 

How,  having  summary  data  by  each  method,  a  comparison  was 
made  region  by  region  for  each  day  to  determine  if  bright  re¬ 
gions  on  the  SOON  scans  corresponded  to  bright  areas  (in 
longitude  only)  on  the  Algonquin  plots.  Since  the  Algonquin 
data  can  only  give  an  indication  of  the  relative  distance  of 
active  regions  from  the  solar  central  meridian,  a  rough  com¬ 
parison  of  longitudes  is  possible  but  not  pinpoint  accuracy  in 
either  latitude  or  longitude.  This  becomes  apparent  when 
looking  at  the  Algonquin  scans.  In  reality  then,  the  latitude 
reading  produced  by  the  SOON  program  is  unnecessary  but  is 
included  for  completeness. 

As  previously  mentioned,  available  data  from  the  SOON 
network  was  limited  to  the  May  to  December  1978  time  frame. 
Of  this,  only  data  from  1  May  through  16  July  1978  could  be 
used  in  the  analysis.  Additional  data  was  unavailable  due  to 
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the  time  and  expense  involved  in  processing  it. 

r~ 


Figure  4.  Algonquin  Solar  Drift  Scan  -  14  Hay  1978 

Table  I 

Sample  SOON  Program  Printout  -  14  May  1978 
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Figure  4  and  Table  I  will  help  in  understanding  the 
methodology  involved.  This  data  set  is  typical  of  all  the 


data  reviewed.  The  solar  drift  curve  is  a  reduction  of  the 
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actual  curve  received  for  14  May  1978.  The  east  and  west 
solar  limbs  are  identified  by  the  tick  marks  labeled  1710  UT 
and  1707  UTf  respectively.  The  USAP  numbered  regions  and 
percentage  of  solar  radii  for  each,  identified  by  the  SOON 
data,  are  also  shown.  CM  is  the  approximate  location  of  the 
central  meridian.  The  summary  table  for  14  May  1978  is 
representative  of  the  output  from  the  SOON  Program.  It  iden¬ 
tifies  the  USAP  numbered  regions  observed  each  day  (REGION), 
their  latitude  and  longitude  in  the  heliocentric  frame  of  the 
sun  (LAT  and  LONG),  total  area  of  all  histograms  of  the  same 
region  (TOTAL  AREA),  number  of  histograms  recorded  for  each 
region  (OBSV),  the  average  area  for  the  region  observed  (AVG 
AREA),  the  average  area  corrected  for  foreshortening  (COR 
AREA),  the  new  latitude  and  longitude  in  the  Algonquin 
reference  frame  (NLAT  and  NLONG) ,  and  the  percent  of  solar 
radius  where  each  region  lies  in  reference  to  the  solar  cen¬ 
tral  meridian  (SOLRAD). 

One  notable  characteristic  is  the  second  highest  spike  on 
the  Algonquin  plot  just  to  the  east  of  the  central  meridian. 
Throughout  all  of  the  scans  that  tracked  this  feature  across 
the  sun,  there  was  never  a  corresponding  SOON  region  that 
could  be  matched  with  it.  This  is  to  be  expected  since  there 
are  some  features  detected  by  one  method  that  may  not  have  a 
match  in  the  other  data.  Based  on  scans  that  came  before  this 
one,  the  various  spikes  on  the  Algonquin  plot  could  be  cor¬ 
related  with  active  regions  identified  by  the  SOON  data.  Note 
that  Region  1099  was  eliminated  immediately  since  its 
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longitude  was  far  beyond  the  65  degree  limit  and  resulted  in 
an  apparent  solar  radius  of  greater  than  90  degrees.  In  com¬ 
paring  the  above  data  sets  it  is  apparent  that  the  SOON  data 
indicates  large  bright  areas  at  73  and  55  percent  of  a  solar 
radius  east  of  the  central  meridian  and  12,  30,  and  70  percent 
of  a  solar  radius  west  of  the  central  meridian.  To  arrive  at 
the  locations  on  the  10.7  cm  radio  drift  scan  plot  for  these 
points,  the  distance  between  the  indicated  east  and  west  limbs 
was  divided  in  half,  and  these  percentages  were  applied  to 
this  figure.  This  procedure  was  considered  reasonable  since 
the  SOON  data  was  unitized  to  solar  radii  by  setting  R^»  1.0 
as  described  in  Appendix  E.  As  theorized,  spikes  correspon¬ 
ding  to  each  of  these  identified  regions  exist  on  the  Algon¬ 
quin  plot. 

This  method  of  rough  comparison  substantiates  the 
existence  of  a  correlation  which  appeared  between  a  majority 
of  the  daily  scans  and  SOON  summary  printouts.  From  the  data 
on  successive  days,  a  given  active  region  could  be  tracked 
across  the  sun  at  roughly  the  correct  solar  rotational 
velocity  for  the  region's  latitude  (Ref  17:31).  The  fact  that 
both  the  Algonquin  scans  and  SOON  data  depicted  this  movement 
further  reinforces  the  existence  of  a  relationship  between  the 
phenomena.  A  point  to  remember  about  the  Algonquin  data  is 
that,  due  to  the  sensitivity  of  the  equipment,  it  must  be  ad¬ 
justed  to  keep  each  day's  reading  from  being  too  great. 
Therefore,  if  a  very  bright  region  is  in  the  scan,  the  gain 
must  be  reduced  to  keep  that  region  within  the  bounds  of  the 
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equipment.  When  the  active  region  rotates  out  of  the  field  of 
view  of  the  instrument  and  the  instrument  gain  is  increased, 
specific  features  may  become  apparent  that  had  been  obscured. 

The  calibration  constants  in  Appendix  B  (Table  II)  make 
all  the  10.7  cm  drift  scan  readings  comparable  and  were  ap¬ 
plied  before  data  analysis.  This  was  done  by  measuring  the* 
peak  in  centimeters  (scale  not  important)  corresponding  to  an 
active  region  identified  by  the  SOON  equipment.  The  values 
obtained  were  then  divided  by  the  calibration  constant  for  the 
day  under  study  to  arrive  at  the  numerical  input  from  the  Al¬ 
gonquin  data  that  would  be  compared  with  its  corresponding 
region  area  as  defined  by  the  SOON  data. 

Da  ta  Analysis 

Once  both  data  sets  were  reduced  to  this  point,  they  were 
in  a  form  for  direct  comparison  to  determine  if  a  relationship 
existed  between  them  with  the  hope  of  finding  a  definite 
mathematical  relationship.  If  such  a  relationship  could  be 
established,  it  could  then  be  used  for  developing  more  ac¬ 
curate  forecasting  models  of  the  10.7  cm  solar  flux.  This 
then,  is  the  major  aspect  of  this  study;  to  determine  if  a 
relationship  exists  between  the  irradiance  of  the  active  re¬ 
gions  identified  by  the  SOON  data  and  the  features  identified 
by  the  Algonquin  plots.  To  attempt  to  establish  this 
relationship,  the  magnitude  of  the  Algonquin  plot  peak  values 
corresponding  to  active  SOON  regions  were  measured  as  men¬ 
tioned  previously.  The  SOON  corrected  areas  made  up  the  other 
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half  of  the  data.  Areas  were  used  since  they  represent  the 
relative  brightness  of  each  region  above  the  quiet  sun  level 
that  corresp.ond  to  active  regions  on  the  Algonquin  plots.  All 
told,  there  were  180  data  points  (individual  features  in  both 
the  Algonquin  and  SOON  data)  available.  These  represented  53 
distinct  USAF  numbered  regions  in  a  two  and  a  half  month 
period.  Of  these,  38  SOON  data  points  and  their  corresponding 
Algonquin  plot  values  had  to  be  discarded  because  they  were 
too  far  into  the  solar  limbs  (greater  than  80  degrees  east  or 
west  of  the  solar  central  meridian).  Even  though  three  of  the 
data  points  had  longitudes  greater  than  65  degrees,  they  were 
left  in  because  they  were  only  slightly  greater  than  65 
degrees . 

Model  Development 

The  final  pairing  resulted  in  142  points  for  input  into 
the  data  analysis  programs  (see  Appendix  F).  Though  more 
would  have  been  desired,  this  sample  should  be  enough  to 
validate  or  refute  the  existence  of  a  relationship  between  the 
SOON  and  Algonquin  data.  Since  there  was  only  one  feature 
from  each  of  the  data  sets,  a  simple  linear  regression  routine 
was  used  to  compare  the  data  (Ref  5:1-22  and  33:320-349,  293- 

300).  The  basic  linear  regression  model  used  was 

Y  *  a  +  $X  +  c  (8) 

where  a  and  ft  are  the  intercept  and  slope,  respectively,  and  e 

2 

is  the  random  error  with  mean  zero  and  variance  o  .  The  model 
further  assumes: 
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1 .  The  sample  is  random, 

2.  The  dependent  variable,  Y  (Algonquin  scan  value),  is 
drawn  from  a  normal  distribution  for  any  given  value  of  X 
(corrected  area  from  SOON), 

3.  The  regression  of  X  and  Y  is  linear, 

4.  All  of  the  Y  values  have  the  same  variance,  and 

5*  The  random  errors,  e^,  are  independent  of  each  other. 
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Figure  5*  Scattergram  of  CALSCN  versus  S00NAREA 
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The  first  attempt  at  a  regression  equation  relating  X  and  Y 
was  apparently  unsuccessful.  The  scattergram  plot  of  the  data 
seen  in  Figure  5  does  not  indicate  any  substantial  relation¬ 
ship  between  SOON  area  and  Algonquin  intensity.  The 
mathematical  formula  determined  gives  even  less  hope  of  a 
relationship: 

CALSCN  =  B  (SOONAREA)  +  A  (9) 

or 

CALSCN  =  .00015  (SOONAREA)  +  1.09429  (10) 

where 

CALSCN  =  calibrated  Algonquin  scan  value,  and 
SOONAREA  =*  corrected  area  from  the  SOON  histograms. 

The  correlation  coefficient,  r  ,  was  equal  to  -59932  with  an 
2 

r  of  *35919  •  The  correlation  coefficient  indicates  how  well 

a  regression  line  characterizes  the  linear  relationship 

2 

between  two  variables  while  r  is  a  measure  of  the  proportion 
of  variance  in  one  variable  explained  by  the  other  and 
measures  the  strength  of  the  relationship  between  the  two 
variables  (Ref  33:279). 

At  this  point,  it  seemed  reasonable  that  a  better 
mathematical  relationship  could  be  found.  Relying  on  tech¬ 
niques  described  by  Neter  and  Vasserman  (Ref  31:131-136),  the 
plot  depicted  in  Figure  5  seemed  to  indicate  that  the  variance 
increases  as  X  (SOONAREA)  increases.  Therefore,  using  their 
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technique  for  cases  in  which  the  standard  deviation  of  the 
error  terms  in  the  regression  equation  is  proportional  to  Xf 
the  following  transformations  were  developed: 

Y'  =  Y/X  (11) 

X'  =  1/X  (12) 

so  that 

Y’  =  a  +  eX'  +  c  (13) 

becomes  the  basic  equation.  In  the  two  analysis  programs, 
Equations  ( 1 1 )  and  ( 1 2 )  become 

NEVSCN  =  CALSCN/SOONAREA  (14) 

and 

NEWAREA  =  1 / S  OONAREA  .  ( 1 5 ) 

When  these  values  were  used  in  the  SPSS  regression  routine  a 
better  correlation  was  established  as  shown  in  the  scattergram 
in  Figure  6.  From  this,  the  resulting  equation  for  predicting 
an  active  region  on  the  Algonquin  scans  with  SOON  data  is: 

NEWSCN  -  B  (NEWAREA)  ♦  A  (16) 

which,  in  the  variables  of  interest,  becomes 

CALSCN/SOONAREA  -  B  (l/SOONAREA)  ♦  A  .  (17) 
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Figure  6.  Scattergram  of  NEVSCN  versus  NEWAREA 


The  correlation  coefficient  for  this  relationship 
giving  an  r  of  .78338  which  indicates  that  the 
equation  defines  a  good  relationship  between  the 
(calibrated  Algonquin  scans)  and  independent  (SOON 
area)  variables.  Equation  (17)  can  be  simplified  to 


is  .88509 
regression 
dependent 
corrected 


CALSCN  -  B  +  A  (SOONAREA) 


(18) 
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resulting  in 


CALSCN  -  .88782838  ♦  .00025237604  (S00NAREA)  ( 1 9) 
using  the  statistics  provided  in  the  SPSS  program. 

Conf idence  Interval 

Equation  (19)  is  good  for  point  estimates  hut  because  the 
deviations  from  the  line  increase  with  large  SOONAREA,  the 
certainty  associated  with  the  point  estimate  becomes  less  for 
large  SOONAREA.  Therefore,  in  order  to  quantify  the  potential 
errors  involved  in  using  this  equation,  confidence  limits  must 
be  established.  The  95  percent  prediction  interval  for  values 
predicted  with  Equation  (19)  is  found  as  follows  (Ref  21:373- 
374): 

V*  -  to/2,n-2/MSE’^1  +  l/n  +  (X1  -  T’^/S^,]  <  YQ/X 
<  Yo/X  ♦  ta/2,h-2/MSE^1  +  l/n  +  (X’  "  (20) 

where 

Yq/X  35  Y  1  *  NEWSCN  value  predicted  using  Equation  (16), 
Yq  «  CALSCN  value  predicted  using  Equation  (19)> 

X  *  SOONAREA  observed, 

Yq  *  The  true  value  of  CALSCN, 

*a/2,n-2  "  *.05/2,140  *  1,97» 

HSE’  *  Mean  square  error  of  the  transformed  data, 

X'  -  Transformed  X  value  for  which  CALSCN  value  is  to  be 
predic  ted , 
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\  -  *./2,n-2/,,SE't'  *  *  <*’  -  X  <  To 

<  Yq  ♦  ta/2,n-2/MSE^1  +  l/n  +  ( x  ’ )  2/sxx .]  *  (21  ) 

with 

MSB'  =  .00019012  , 

X*  =  .0007445  ,  and 

Sxx'  =  141  C*000^?06)2  =  .0001221083  , 
so  the  confidence  interval  becomes 

Yo-1 . 97^.0001 901 2[l +1 /I 42+(X’ -. 0007445 )2/ • 0001 221083]  X  <  Yq 

<  Y  *1  . 97/. 0001 901 2[ 1+1 /I  42+ (X' -.0007445) */. 0001 221 083]  X  .  (22) 

Equation  (22)  verifies  that  the  width  of  the  interval  in¬ 
creases  as  X  increases  which  corresponds  to  what  is  observed 
in  the  original  data. 

Equations  (19)  and  (22)  are  consistent  with  the  five  as¬ 
sumptions  expected  of  a  linear  regression  model.  The  results 
indicate  the  sample  was  random,  the  dependent  variable,  Y  , 
came  from  a  normal  distribution  for  any  given  value  of  X,  and, 
after  the  transformation,  the  regression  was  linear,  the 
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variation  of  Y  was  constant,  and  the  random  errors  were  in¬ 
dependent* 

Other  Relationships 

After  reviewing  the  residual  analysis  performed  on  the 
data  by  the  analysis  program  (See  Appendix  P),  it  seemed  that 
at  least  one  outlying  data  point  may  have  been  unduly  influen¬ 
cing  the  regression  equation*  All  together,  ,  there  were  six 
residuals  greater  than  two  standard  deviations.  These 
represented  about  five  percent  of  the  total  and  were  expected, 
however,  one  point  (CALSCN  =  3*04,  SOONAREA  *  668.29)  had  a 
residual  twice  as  great  as  the  next  largest  residual.  After 
rechecking  the  data  reduction  for  this  point  and  finding 
nothing  wrong,  it  was  decided  to  explore  the  effects  of  this 
point  on  the  regression  equation.  The  point  was  eliminated 
from  the  data  and  the  programs  were  rerun.  Equations  ( 1  1  ) 
through  ( 1 8 )  still  apply  only  now  the  resulting  equation  is: 

CALSCN  =  .869333787  +  .00024499284  (SOONAREA)  (23) 

p 

which  resulted  in  an  r  of  *91565  and  r  of  .83842  on  the 
transformed  equation.  The  confidence  limits  for  values 
predicted  with  this  equation  are: 

Y^  -  1  .97/.0001 27?6[ 1 +1 /I  4 1 + ( X ' - . 0007391  ) */ . 0001 21 5552]  X  <  Y 

<.  Yo+1 •97/.0001 2736[l +1/141 +(X’-. 0007391 )2/.0001 21 5552]  X  (24) 

which  has  a  derivation  similar  to  Equation  (22).  Even  though 
the  linear  association  between  the  dependent  and  independent 
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variables,  as  measured  by  the  correlation  coefficient,  is  bet¬ 
ter  than  the  original  formulation,  there  is  no  valid  reason  to 
exclude  the  data  point  in  question  other  than  the  fact  that  it 
is  not  typical  of  the  rest  of  the  data*  For  this  reason,  the 
model  defined  by  Equations  (19)  and  (22)  is  the  best  relation¬ 
ship  found  between  the  Algonquin  scans  and  the  SOON  data* 
Note  also  that  even  with  the  point  removed,  Equation  (23)  does 
not  differ  appreciably  from  Equation  (19)* 

To  further  explore  the  correspondence  between  the  Algon¬ 
quin  scan  values  and  the  SOON  data,  a  second  set  of  data  was 
developed  from  the  original  histogram  records.  Now,  however, 
absolute  brightness  was  calculated  for  each  •  90N  region  by 
multiplying  each  area  by  its  bin  value  (relative  brightness). 
Since  each  bin  is  linearly  brighter  than  the  preceeding  one, 
this  should  give  a  good  indication  of  how  active  the  region 

is.  This  means  that  the  larger  a  bin  area  is,  the  greater  its 

\ 

contribution  to  the  overall  brightness  of  the  region  which,  in 
turn,  should  be  reflected  by  a  larger  spike  in  the  correspon¬ 
ding  Algonquin  scan.  As  expected,  similar  results  were  ob¬ 
tained  with  this  method  as  in  the  first  analysis  and  will  not 
be  pursued  further  here. 
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IV.  Conclusions  and  Recommends tions 


Conclusions 

From  this  analysis  it  is  apparent  that  a  relationship 
does  exist  between  Algonquin  drift  scans  of  the  sun  and  SOON 
histogram  records  of  active  regions  on  th§  sun.  The  area 
above  the  quiet  sun  level  calculated  from  the  SOON  data  can  be 
used  to  predict  the  location  and  intensity  of  solar  bursts  in 
the  10.7  cm  solar  radio  range  as  gleaned  from  the  Algonquin 
scans  of  the  sun.  With  further  analysis,  as  recommended 
below,  the  SOON  histograms  could  become  a  valuable  tool  in 
aiding  the  understanding  of  the  relationship  between  solar 
active  region  area  and  the  10.7  cm  solar  flux. 

Recommendations 

In  order  to  improve  this  analysis  several  important 
modifications  to  this  experiment  could  be  done  to  refine  the 
resul ts • 

1 .  The  ephemeris  data  used  in  the  analysis  could  be 
modified  to  account  for  the  fact  that  the  tabulated  values  are 
for  0000Z  while  the  Algonquin  solar  scans  and  SOON  records 
were  taken  at  1700Z  each  day.  Though  these  changes  are  small, 
at  certain  times  of  the  year  the  angles  do  vary  rapidly  and 
could  affect  the  accuracy  of  the  data. 

2.  To  use  as  much  data  as  possible,  corrections  could  be 
applied  to  those  SOON  data  records  with  a  wavelength  offset 
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other  than  zero.  Though  few  records  have  a  wavelength  offset 
other  than  zero,  those  that  do  could  be  significant  and  should 
not  be  wasted. 

3.  Further  discussions  with  the  SOON  engineers  and 
physicists  might  uncover  better  ways  to  categorize  and  quan¬ 
tify  the  histograms.  In  this  study,  the  quiet  sun  value  was 
used  as  a  benchmark  to  determine  what  brightness  levels  (bin 
numbers)  should  be  summed  to  get  a  rough  idea  of  a  region's 
brightness  above  the  background  value,  thus  making  an 
inference  as  to  the  region's  contribution  to  hydrogen-alpha 
wavelength,  and,  after  regression,  10.7  cm  flux.  A  better 
benchmark  might  be  found  such  as  using  the  flare  threshold 
rather  than  brightness  above  quiet  sun. 

4-  Data  used  in  any  follow  on  study  must  be  standardized. 
The  only  data  available  for  this  study  was  in  the  May-December 
1978  time  frame.  The  SOON  system  was  very  new  at  the  time  and 
was  going  through  considerable  growing  pains.  This  is  reflec¬ 
ted  in  the  fact  that  data  was  not  all  standarized  to  the  same 
quiet  sun  level  (some  was  "standarized"  to  bin  10  while  some 
was  "standarized"  to  bin  12).  If  the  same  time  frame  is  used, 
the  data  must  all  be  standarized  to  the  same  level. 

5*  Though  the  Algonquin  scans  are  supposed  to  be  taken  at 
1700Z,  they  often  slip  one  way  or  the  other  by  as  much  as  20 
minutes.  Since  the  SOON  data  was  only  available  at  1700Z,  if 
this  was  the  case,  the  data  was  compared  as  if  both  readings 
were  taken  at  1700Z.  In  the  future,  the  times  of  both  data 
sets  should  be  matched  exactly. 
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6*  A  study  such  as  this  might  be  expanded  to  try  and  find 
a  relationship  between  SOON  data  and  other  wavelengths*  Such 
studies  would  involve  other  SOON  sites  besides  just  Holloman, 
New  Mexico. 

7*  In  the  future,  data  must  be  analyzed  for  a  longer 
period.  Though  data  from  May  through  December  1978  was 
available,  not  all  of  it  could  be  used.  Future  studies  should 
expand  to  cover  many  months  (6  or  more)  and  even  years  rather 
than  just  2  1/2  months,  especially  when  many  days  are  lost  to 
weather  and  equipment  problems.  Additionally,  as  many  regions 
as  possible  should  be  included  each  day.  Finally,  studies 
should  explore  any  relationship  between  the  solar  cycle  and 
these  results,  looking  specifically  at  periods  during  solar 
maximum  and  solar  minimum. 

8.  A  more  complex  method  should  be  investigated  for 

determining  the  intensity  of  the  Algonquin  data  than  measuring 

\ 

the  peak  values.  One  possibility  is  to  take  an  average  value 
for  each  active  longitude  (as  identified  by  a  spike).  Another 
would  be  to  integrate  the  area  under  the  curve. 

9.  Analysis  should  be  done  to  perfect  methods  to  predict 
the  10.7  cm  solar  flux  using  the  Algonquin  solar  scans.  Once 
this  is  done,  based  on  the  present  study,  SOON  data  could  be 
used  to  predict  the  10.7  cm  solar  flux. 
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APPENDIX  A 


Space  Environment  Services  Center  Forecast 


The  following  example  is  the  primary  forecast  issued  on 

19  July  1974.  Note  the  10.7  cm  forecast  and  90-day  mean  in 

Section  IV  (Ref  28:49)* 

NNNN 

KBOU 

HFUS  BOU  192200 

FROM  SPACE  ENVIRONMENT  SERVICES  CENTER  BOULJDER  COLO 
SDF  NUMBER  200 

JOINT  AFGWC/SESC  PRIMARY  REPORT  OF  SOLAR  AND  GEOPHYSICAL  ACTIVITY  ISSUED 
2200Z  19  JULY  1974 

IA.  SOLAR  ACTIVITY  HAS  EEEN  VERY  LOW  WITH  THREE  NON-ENERGETIC  S'JBFLARES 
REPORTED  DURING  THE  PAST  24  HOURS.  REGIONS  438  (S08W72)  AND  443  (SL2EC6) 
HAVE  BEEN  STABLE.  MINOR  INTENSITY  FLUCTUATIONS  HAVE  OCCURRED  IN  REGIONS 
442  (SL0W34)  AND  445  (SC4W37).  AN  ACTIVE  PROMINENCE  (SL3E90)  PARTIALLY 
ERUPTED,  AND  THEN  TOTALLY  DISSIPATED,  BETWEEN  0715-0340Z.  NO  OTHER  ACTI¬ 
VITY  HAS  BEEN  REPORTED  TO  SUBSTANTIATE  AN  ACTIVE  RETURN  OF  OLD  REGION  435. 

IB.  SOUR  ACTIVITY  IS  EXPECTED  TO  REMAIN  LOW. 

II.  THE  GEOMAGNETIC  FIELD  HAS  BEEN  QUIET.  IT  IS  EXPECTED  TO  BE  QUIET  TO 
UNSETTLED. 

III.  EVENT  PROBABILITIES  20  -  22  JULY 
CUSS  M  03/03/03 

CUSS  X  01/01/01 
PROTON  01/01/01 
PCAF  GREEN 

IV.  OTTAWA  10.7  CM  FLUX 
OBSERVED  19  JULY  84 
PREDICTED  20  -  22  JULY  83/85/85 
90-DAY  MEAN  19  JULY  90 

V.  GEOMAGNETIC  A  INDICES 
OBSERVED  FREDERICKSBURG  18  JULY  06 
ESTIMATED  AFR/AP  19  JULY  04/06 

PREDICTED  AFR/AP  20  -  22  JULY  07/08  -  09/10  -  11/12 

SOLTERWARN 

SPAN 

BT 
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APPENDIX  B 


Algonquin  Videometer  Scan  and  Calibration  Levels 

The  following  videometer  scan  (Figure  7)  is  typical  of 
the  data  received  from  the  Algonquin  Radio  Observatory,  These 
same  scans  are  reproduced  in  Solar  Geophysical  Da  ta  (May  1978- 
July  1978)  in  a  much  smaller  format.  The  east  and  west  limbs 
of  the  solar  disk  are  annotated  on  the  graph  and  the  point 
halfway  between  these  points  gives  an  approximate  indication 
of  the  Central  Meridian  location.  Based  on  this  data,  rough 
correlations  with  longitude  could  be  made  with  active  regions 
identified  by  the  SOON  data.  The  spikes  on  this  graph 
represent  longitudes  with  higher  than  normal  activity.  The 
horizontal  baseline  is  for  reference  only  and  had  to  be  ad¬ 
justed  slightly  to  standardize  all  of  the  plots.  Adjustments 
took  the  form  of  adding  or  subtracting  a  small  correction  from 
the  plot  readings  if  the  scan  did  not  start  precisely  on  the 
baseline.  These  readings  were  then  further  standardized  by 
applying  the  daily  calibration  constants  listed  in  Table  2. 
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Figure  7.  Algonquin  Solar  Drift  Scan  -  1  May  1978 
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Table  II 


Calibration  Levels 


May 

June 

July 

Aug. 

Sept . 

Oct . 

Nov. 

Dec . 

1 

3.15 

6.16 

4.99 

5.71 

3.09 

5.59 

5.45 

5.05 

2 

3.04 

5.93 

5.53 

5.41 

2.99 

5.87 

4.99 

5.23 

3 

3.18 

6.37 

5.79 

5.60 

3.04 

5.53 

4.71 

6.23 

4 

3.19 

6.19 

5.71 

5.38 

3.11 

5.27 

4.26 

5.51 

5 

3.03 

6.37 

5.87 

5.74 

3.01 

5.56 

5.33 

4.67 

6 

3.01 

5.86 

5.60 

5.77 

2.87 

5.87 

5.45 

5.19 

7 

2.81 

6.05 

5.56 

5.81 

2.77 

5.56 

4.73 

5.44 

8 

2.76 

5.92 

5.84 

5.62 

2.94 

5.61 

4.38 

5.89 

9 

5.84 

NO  DATA 

5.70 

5.96 

2.82 

2.66 

4.30 

2.13 

10 

5.86 

4.73 

2.83 

6.02 

2.87 

2.53 

4.47 

2.13 

11 

5.43 

6.02 

2.68 

5.72 

5.73 

2.64 

4.45 

1.37 

12 

6.00 

5.98 

2.80 

5.72 

5.68 

4.39 

4.38 

2.71 

13 

6.10 

5.44 

2.74 

5.66 

5.81 

4.90 

5.05 

2.40 

14 

6.29 

5.92 

2.79 

5.51 

5.87 

5.26 

4.04 

2.50 

IS 

5.80 

6.05 

2.68 

5.27 

5.85 

5.13 

4.21 

2.20 

16 

6.06 

5.99 

2.77 

5.08 

6.19 

5.51 

4.60 

2.32 

17 

6.21 

5.81 

2.77 

5.47 

6.08 

5.45 

4.80 

4.98 

18 

5.84 

5.97 

2.70 

5.61 

6.10 

5.33 

4.47 

5.20 

19 

5.66 

5.95 

2.28 

5.75 

3.00 

4.93 

4.20 

5.57 

20 

6.38 

5.42 

5.57 

5.18 

2.97 

5.38 

5.33 

5.25 

21 

6.25 

6.03 

5.28 

5.74 

2.80 

5.35 

5.94 

4.54 

22 

NO  OAT A 

5.52 

5.61 

5.64 

6.21 

5.30 

5.05 

5.31 

23 

6.27 

5.67 

5.66 

5.79 

2.93 

5.19 

4.95 

5.07 

24 

2.78 

5.74 

6.19 

5.61 

2.88 

5.24 

5.65 

5.21 

2S 

5.83 

6.73 

5.97 

5.48 

2.86 

4.56 

4.69 

5.26 

26 

6.19 

2.79 

5.06 

5.75 

3.32 

4.48 

NO  OATA 

2.46 

27 

6.15 

6.19 

5.90 

5.91 

2.86 

4.70 

.  4.83 

2.53 

28 

6.32 

5.76 

5.67 

6.15 

2.71 

4.12 

4.98 

2.47 

29 

5.99 

5.66 

5.65 

6.00 

2.81 

4.96 

4.97 

4.85 

30 

6.15 

5.54 

6.26 

2.85 

5.63 

5.14 

5.07 

5.53 

31 

6.59 

6.03 

6.00 

4.36 

5.28 
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APPENDIX  C 


SOON  Card  Forma t 

The  SOON  data  was  received  on  punched  cards.  A  region 
identifier  record  format  is  shown  first.  These  cards  are  used 
by  the  SOON  scientists  to  identify  regions  to  be  observed  by 
the  SOON  telescopes  and  give  region  identifying  information 
that  is  inherent  in  each  of  the  histograms  which  follow  for 
that  region.  Following  this  is  the  format  for  the  SOON 
histogram  data  sets.  The  basic  card  layout  is  specified  and 
then  an  example  card  image  is  presented.  The  SOON  program 
outlined  in  Appendix  D  follows  this  format  to  reduce  the  data. 


FORMAT  FOR  REGION  IDENTIFIER  RECORDS 


COLUMN  PARAMETER  DESCRIPTION 


1 

ICONT 

CONTINUATION  CODE 

2 

NOBS 

OBSERVATION  CODE 

3-6 

IDENT 

RGN  §  or  RGDF 

7-8 

MTIME(l) 

YEAR 

9-10 

MTIME(2) 

MONTH 

11-12 

MTIME(3) 

DAY 

13-14 

MTIME(4) 

HOUR 

15-16 

MTIME(5) 

MINUTE 

17-18 

MTIME(6) 

SECOND 

19-22 

NIDRGN 

NEW  REGION  IDENTIFIER 

23-24 

IRCHG 

CHANGE 

25-28 

REGID 

REGION  IDENTIFIER 

29-33 

IRBITS 

REGION  CONTROL  BITS 

34-38 

SLONG 

REGION  CENTER  LONGITUDE  (POS=WEST,NEG=EAST) 

39-43 

SLAT 

REGION  CENTER  LATITUDE  (POS=NORTH,NEG=SOUTH) 

44-53 

RV 

BOX  CENTER  RADIUS  VECTOR 

54-63 

SINPA 

SINE  OF  THE  P  ANGLE  FOR  BOX  CENTER 

64-73 

COSPA 

COSINE  OF  THE  P  ANGLE  FOR  BOX  CENTER 

74-83 

DLTLG 

SOLAR  ROTATION  RATE  AT  THIS  LATITUDE  (RAD.) 

84-88 

SIZEEW 

E-W  SIZE  IN  VIDOMETER  UNITS 

89-93 

SIZENS 

N-S  SIZE  IN  VIDOMETER  UNITS 

94-95 

THRSHS 

THRESHOLDS 

96-97 

. .  THRSHQ 

98-99 

THRSHP 

100-101 

THRSHF 

102-106 

AVBIN 

AVERAGE  BRIGHTNESS 

Figure  8.  Format  for  Region  Identifier  Records 
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FORMAT  FOR  HISTOGRAM  RECORDS 


*****  CARD  IMAGE  RECORDS  ***** 

Card  10:  (1st  card  of  sequence) 

OTMN  PARAMETER  DESCRIPTION 


1 

ICONT 

CONTINUATION  CODE 

2 

NOBS 

OBSERVATORY  CODE 

3-6 

IDENT 

REGION  IDENTIFIER 

7-8 

MT IME ( 1 ] 

)  YEAR 

9-10 

MTIME(2; 

)  MONTH 

11-12 

MTIME(3, 

)  DAY 

13-14 

MT IME (4 ; 

)  HOUR 

15-16 

MTIME(5’ 

)  MINUTE 

17-18 

MTIME(6] 

)  SECOND 

19-23 

SLAT 

LATITUDE  OF  REGION  CENTER  (RADIANS) 

24-28 

SLONG 

LONGITUDE  OF  REGION  CENTER  (RADIANS) 

29-32 

DLTWLJ 

WAVELENGTH  OFFSET  (from  H-alpha  center) 

33-36 

setup(i; 

>  STANDARD  OPTICAL  TV  SETUPS 

37-40 

$ETUP(3 

) 

41-44 

SETUP ( 5  * 

> 

45-49 

MODE  2 

STATISTICAL  MODE  OF  SECOND  READ 

50-55 

CAL 

CALIBRATION  CONSTANT 

56-60 

JNF 

NUMBER  OF  FIELDS  INTEGRATED 

61 

BLANK 

62-78 

SORTKEY 

DUPL  OF  COL  2-18  USED  FOR  SORTING  ONLY.. 

79-80 

CARD  NO, 

.  VALUE  *  10 

Figure  9*  Format  for  Histogram  Records 


\  f 

V 

l.  : 


SOON  HISTOGRAM  FORMAT  -  Card  Image 
*****  DETAIL  OF  PARAMETERS  ***** 


COLUMN  definition: 


1  ■  ■■ 

2 

3 

4 

5 

6  7  8 

1  0 

0 

0 

0 

0 

0  0  0 

2H10927805031 70532 

336 

497 

0  642 

1  210 

8857  45722 

10  H10927805031 7053210 

0 

0 

0 

94 

894  16013106751  H1092780503170S3220 

216157397469112164 

51146 

21149 

6161 

2523 

1440 

926 

683  H10927805031 7053230 

492  452  423 

347 

285 

304 

357 

244 

59 

44  H10927305031 7053240 

1 

19 

4 

0 

0 

0 

0 

0  H10927305031 7053 250 

0  0  0 

0 

0 

0 

0 

0 

0 

0  H10927805031 7053260 

0  0  0 

0 

0 

0 

0 

0 

0 

0  H10927805031 7053270 

0  0  0 

0 

0 

0 

0 

0 

0 

0  H10927805031 7053280 

***  DECOOE  *** 


Card  10:  (see  col  79-80) 

2H1092  -  Double  size  Histogram 
Holloman  AFB  SOON  site 
USAF  region  no.  1092 
780503  -  May  3,  1978 
170532  -  17hr  Q5m1n  32sec  UT 
336  -  .336  radians  =  solar  lat 
=  19.2  north 
497  -  .497  radians  *  solar  CMD 
>  28.5  west 

0  -  H-Alpha  center  line  (6562,8) 
642  -  Setup(l)  equlment  conflg. 

1  -  2  * 

210  3 

8857  -  8.857  or  approx  bln  9*Q  Sun 
45722  -  Calibration  const. 

10  -  JNF 


Card  20: 

- r- 

0  - 
0  - 
94  - 
894  - 
16013  - 
106751  - 
Card  30: 
715157“- 


683 

Card  40: 


B1n(l)  area  MSD  integrated 
over  entire  box 
B1n(2)  ■ 

B1n(3)  MSD  »  Millionths  of 
B1n(4)  solar  disk 

81n(5)  uncorrected 

Bin(6)  • 

Bin(7) 


Bln (8)  M 
(17)  - 


B1n(18)  ■ 

fcr 


******  NOTE:  BRIGHTNESS  BINS  ABOVE  27  ARE  ZERO  FOR'  QUIET-TIME  HISTOGRAMS  ****** 


Card  SO: 

19  -  81 n(28) 
0  -  B1n(34) 

Card  60: 

-  0  -  B1n(3S) 

0  -  B1n(44) 


Card  70: 

iT-  B1n(45) 

0  -  BIn(54> 

Card  80: 

IT-  B1n(55) 

0  .  B1n(64) 


***  NOTE:  CAROS  50-80  ARE  ONLY  USEO  IF  BRIGHTNESS  BIN(27)  IS  EXCEEDED  *** 


Figure  10.  SOON  Histogram  Format  -  Card  Image 
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APPENDIX  D 


SOON  Program 

The  following  program  reduced  the  SOON  histogram  inputs 
into  a  usable  and  comparable  form.  The  data  was  received  on 
punched  cards  in  the  format  depicted  in  Appendix  C.  The  hun¬ 
dreds  of  data  records  involved  were  run  through  the  program 
and  outputted  in  the  tabular  format  described  in  the  text. 
Comments  are  included  in  the  program  to  document  what  is  oc¬ 
curring  at  each  point.  Many  of  the  operations  performed  may 
seem  redundant  but  they  were  left  in  this  format  so  each  could 
be  readily  identified  and  as  an  aid  to  clarity.  Appendix  E 
contains  a  detailed  explanation  of  the  coordinate  transforma¬ 
tion  used  to  get  the  data  into  a  common  reference  frame. 
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1 00-CAP, CM70000, Ti 00 , 10120.  T820 740, PUZ ,4172 

110-FTN5. 

120-LCO. 

1 30-*E0R 

140-  PROGRAM  SOON 

150-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


160-C  THIS  PROGRAM  REDUCES  THE  SEMI-PROCESSED  SOON  DATA  INTO  C 
170-C  DAILY  SUMMARIES  OF  SOLAR  ACTIVE  REGIONS.  THESE  REGIONS  C 
180-C  WILL  BE  IDENTIFIED  BY  USAF  REGION  NUMBER,  SOLAR  LATITUDE  C 
190-C  AND  LONGITUDE  BASED  ON  THE  HELIOCENTRIC  CENTER  OF  THE  SUN  C 
200-C  FOUND  IN  THE  AMERICAN  EP11EMERIS  AND  NAUTICAL  ALM  ANAC- 1  9  7  8  ,  C 
210-C  BRIGHT  REGION  AREA,  AND  REVISED  LATITUDE  AND  LONGITUDE  IN  C 
220-C  THE  REFERENCE  FRAME  OF  THE  ALGONQUIN  DATA.  C 
230-C  CAPT  CRAIG  A.  PUZ  AFIT/ENA  GSO  82D  23  OCTOBER  1982  C 


240-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

250-C 

260-  INTEGER  YRMODA , T IME , AREA , I , J , M , N , DATA(  64 ) , FLAG, FOUND, 

270-  +LPCNTR, YESTDY , EPH DAY , DAY , QB I N , QS UN , LA.LN, NOREG , WAVOF F , CALI 

280-  REAL  LAT , LON, PI, LOS, Z, 

290-  +ALFA(245,4) .STORE (14, 9) , RAD , THETA , BO , LO , A , B , C , P , CAL 

300-  CHARACTER  RECI0N*4 , RECONS ( l 4) *4 

310-C 

320-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
330-C  INITIALIZE  ALL  VARIABLES  AND  ARRAYS.  EPH  DAY  SPECIFIES  THE  C 


340-C  NUMBER  OF  DAYS  THE  PROGRAM  WILL  RUN.  IT  MUST  BE  UPDATED  C 
350-C  PRIOR  TO  EACH  RUN.  YESTDY  IS  THE  DATE  OF  THE  FIRST  DATA  C 
360-C  SET  (YEAR/MONTH/DAY  TWO  DIGITS  EACH).  LAT  AND  LON  ARE  C 
370-C  THE  LATITUDE  AND  LONGITUDE  OF  THE  FIRST  REGION  SCANNED  ON  C 
380-C  ON  THE  FIRST  DAY.  C 


390-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


400-C 

410- 

AREA-0 

420- 

LPCNTR-0 

430- 

YEST  DY-7  80  50 1 

440- 

LAT— 481 

450- 

LON-1266 

460- 

PI-ACOS(-l.O) 

470- 

RAD-2  *P I 

480- 

NOREC-O 

490- 

500-C 

EPHDAY-9 

510- 

DO  3  1-1,64 

520- 

DATA( I ) -0 

530- 

540-C 

3 

CONTINUE 

550- 

9—4 

—4 

1 

t-i 

o 

CM 

O 

o 

560- 

REGONS ( I ) -  * 

570- 

580-C 

20 

CONTINUE 

590- 

DO  7  1-1,14 

600- 

DO  9  J-1,9 

610- 

STORE(I, J)-0. 

620- 

9 

CONTINUE 

630- 

640-C 

7 

CONTINUE 

650- 

DO  6  N-l , EPHDAY 

660- 

DO  4  M-l , 4 

670- 

ALF A( N , M) -0 . 0 

680- 

4 

CONTINUE 

690- 

700-C 

6 

CONTINUE 
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710-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
720-C  READ  E PUEMERI S  DATA  IN  TO  AN  ARRAY.  P-POSITION  ANCLE  OF  C 
730-C  THE  NORTHERN  EXTREMITY  OF  THE  AXIS  OF  ROTATION  OF  THE  SUN  C 
740-C  MEASURED  EASTWARD  FROM  THE  NORTH  POINT  OF  THE  DISK.  BO  AND  C 
750-C  LO-HEHO  GRAPHIC  LATITUOE  AND  LONCITUDE,  RESPECTIVELY,  OF  C 
760-C  THE  CENTER  POINT  OF  THE  SOLAR  DISK.  C 
770-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


780-C 

790- 

DO  15  N»1 , EPHDAY 

800- 

READ( *  ,  5)  (ALFA(N,H) ,M-1 ,4) 

810- 

5 

FORMAT(1X,F7.0,1X,F6.2,2XIFS.2,2X,F6.2) 

820- 

830-C 

15 

CONTINUE 

840- 

DO  12  N-l, EPHDAY 

850- 

WRITE(*  ,8)  (ALFA(N.M) ,M-1 ,*) 

860- 

12 

CONTINUE 

870- 

880-C 

8 

F0RMAT(1X,F7.0,1X,F6.2,2X,F5.2,2X,F6.2) 

890-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


900*C  READ  FIRST  CARD  OF  NEW  DATA  SET  (HISTOGRAM).  FLAG-1  OR  2  C 
910-C  SPECIFYING  SINGLE  OR  DOUBLE  LENGTH  HISTOGRAM.  YRMODA-YEAR ,  C 
920-C  ffONTH ,  AND  DAY  OF  HISTOGRAM.  TIME-HOUR,  MINUTE,  AND  SECOND  C 
930-C  OF  READING  (NOT  USED).  LA , LN-REG ION  LATITUDE  AND  LONGITUDE  C 
940-C  IN  RADIANS  (DIVIDED  BY  1000  FOR  ACTUAL  READING).  WAVOFF-  C 
950-C  WAVE  LE  NOTH  OFFSET  (PROXIMITY  OF  READING  TO  HYDROGEN-ALPHA  C 
960-C  CENTERLINE).  HISTOGRAM  IS  DISCARDED  IF  THIS  IS  OTHER  THAN  C 
970-C  ZERO.  QSUN-QU IET  SUN  BIN  VALUE  ON  HISTOGRAM.  AFTER  C 

980-C  ROUNDING  TO  THE  NEAREST  INTEGER,  ALL  BINS  ABOVE  QSUN+2  ARE  C 
990-C  ADDED  AND  REPRESENT  THE  RELATIVE  BRIGHTNESS  OF  THE  REGION.  C 
1000-C  CALI  IS  A  CALIBRATION  CONSTANT.  C 


1010-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

1020*0 

1030-  1  READ( *  , 10, END-300 ) FLAG , REGION, Y RMO DA, TIME, LA, LN, WAVOFF, 

1040-  +QSUN , CALI 

1050-  10  FORMAT (II, IX , A4 , 216 , 21 5 , 14 , 1 2X , I 5 , 16) 

1060-C 

1070-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1080-C  CONVERT  LATITUDE  AND  LONGITUDE  TO  REAL  VALUES  IN  RADIANS.  C 
1090-C  CONVERT  CALIBRATION  CONSTANT  INTO  REAL  VALUE.  C 

1100-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

lilO-C 

1120-  LAT-REAL( LA) /i 000.0 

1130-  LON-REAL ( LN)/1000.0 

1140-  CAL-REAL(CALi )/l .0E6 

1150-C 

1160-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1170-C  COMPARE  DATE  ON  LATEST  CARD  WITH  PREVIOUS  CARDS  DATE.  IF  C 

1180-C  IT  IS  THE  SAME,  CONTINUE  ADDING  TO  PRESENT  DAY.  OTHERWISE,  C 

1190-C  IT  IS  A  NEW  DAY  AND  THE  PREVIOUS  DAY’S  DATA  MUST  BE  C 

1200-C  SUMMARIZED.  C 

1210-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1220-0 

1230-  2  IF  (YESTDY.EQ.YRMODA)  THEN 

L240-C 

I  2  50-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1260-C  DISTINGUISH  BETWEEN  DOUBLE  AND  SINGLE  SIZE  HISTOCRAMS.  C 

1270-C  READ  APPROPRIATE  NUMBER  OF  DATA  CARDS.  C 

1280-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

1290-C 
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1300-  IF  (FLAC.EQ.2)  THEM 

1310-  LPCNTR-6  4 

1320-  READ<*,30)  ( DATA( I ) # 1-1 , 7 ) 

1330-  READ ( * , 32 )  (DATA( I) ,1-8, 17) 

1340-  READ( *  ,  32  )  < DATA( I ) , 1-18 t 27  ) 

1350-  READ( *  t  30)  ( DATA( I ) , 1-28 B 34) 

1360-  READ( * , 32  )  ( DATA( I ) # 1-3 5 , 44  ) 

1370-  READ( * , 32 )  ( DATA( I ) , I -4 5 , 54  ) 

1380-  READ( * ,  32  )  ( DATA( I ) # 1-5 5 # 64  ) 

1390-  30  F0RMAT( 1 8X  t  7 1 6  ) 

U00-  32  FORMAT  (1016) 

1410-  ELSE 

1420-  LPCNTR-2  7 

1430-  READ ( * . 30)  (DATA( I ) , 1-1  *  7 ) 

1440-  READ( *  ,  32  )  (DATA(I) 1 1  -8  , 17) 

1450-  READ<*,32)  (DATA(I) ,1-18,27) 

1460-  END  IF 

1470-C 

1480-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


1490-C  CHECK  WAVELENGTH  OFFSET.  IF  IT  EQUALS  2ER0  DETERMINE  THE  C 
1500-C  QUIET  SUN  BIN  LEVEL  (QBIN)  AND  ADD  THE  AREAS  OF  ALL  BINS  C 
1510-C  ABOVE  THIS  VALUE  (QBIN+2  TO  THE  NUMBER  OF  BINS  READ-E ITHER  C 
1520-C  27  OR  64).  PLAGE  GENERALLY  BEGINS  TWO  BINS  ABOVE  QBIN.  C 


t 530-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1540-C 

155^-  IF  (WAVOFF.EO.O)  THEN 

1560-  QBIN-NINT( REAL (QSUN)/1000.0) 

1570-  QBIN-QBIN+2 

1580-C 

1590-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1600-C  SUM  THE  AREA  OF  ALL  BINS  ABOVE  THE  QUIET  SUN  LEVEL.  C 

1610-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1620-C 

1630-  DO  60  J-QBIN , LPCNTR 

1640-  AREA-DATA(J )+ARE A 

1650-  60  CONTINUE 

1660-C 

1670-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1680-C  REJECT  THE  ENTIRE  HISTOGRAM  SINCE  WAVELENGTH  OFFSET  IS  C 

1690-C  OTUER  THAN  ZERO.  C 

1700-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1710-C 

1720-  ELSE 

1730-  GO  TO  1 

1740-  END  IF 

1750-C 

1760-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1770-C  INITIALIZE  AND  INCREMENT  COUNTER  FOR  DATA  STORAGE  ARRAY.  C 
1780-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1790-C 

1800-  1-0 

1810-  100  I-I+l 

1820-C 


1830-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1840-C  FOR  EACH  REGION  OF  A  GIVEN  DAY,  ASSIGN  THE  FOLLOWING  VALUES  C 


1850-C  TO  EACH  COLUMN:  C 
1860-C  COL  O-USAF  REGION  NUMBER  C 
1870-C  COL  1-LATITUDE  OF  REGION  VIEWED,  C 
1880-C  COL  2-LONGITUDE  OF  RECION  VIEWED,  C 
1890-C  COL  3-TOTAL  AREA  OF  BRIGHT  REGIONS  ABOVE  QUIET  SUN  LEVEL  C 
1900-C  ( QBIN)  IN  MILLIONTHS  OF  TUE  SOLAR  DISK,  C 
1910-C  COL  4-COUNTER  FOR  NUMBER  OF  TIMES  EACH  REGION  IS  SUPPLE-  C 
1920-C  MENTED  BY  ANOTHER  HISTOGRAM,  C 
1930-C  COL  5-AVERAGE  AREA  FOR  EACH  USAF  REGION  OBSERVED,  C 
1940-C  COL  9-PERCENT  OF  SOLAR  RADIUS  (FROM  E-W  CENTRAL  MERIDIAN)  C 
1950-C  WHERE  ACTIVE  REGION  LIES.  C 


1960-C  NOREG-COUNTER  FOR  NUMBER  OF  DIFFERENT  REGIONS  READ  PER  DAY.  C 
1970-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
1980-C 

1990-  IF  (ST0RE(I,3).EQ.0.0)  THEN 

2000-  RECONS ( I )-REG ION 

2010-  STORE ( 1,1 ) -LAT 

2020-  STORE (I,2)-L0N 

2030-  STORE (I , 3 ) -REAL< ARE A ) *CAL 

2040-  STORE ( I , 4 ) *1 

2050-  FOUND- l 

2060-  NOREG-NOREG+i 

2070-  ELSE 

2080-C 

2090-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2100-C  A  CHECK  FOLLOWS  TO  DETERMINE  IF  THE  LATEST  REGION  IS  C 

2110-C  ONE  PREVIOUSLY  MEASURED  OR  A  NEW  ONE.  IF  IT  IS  A  NEW  C 

212 0-C  REGION,  ANOTHER  ROW  IS  STARTED  IN  THE  STORAGE  ARRAY.  C 

2130-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2140-C 

2150-  IF  (RECION. EQ . RECONS ( I ) )  THEN 

2160-C 

2170-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2180-C  IF  THE  CURRENT  REGION  DOES  NOT  DIFFER  FROM  ONE  PREVIOUSLY  C 
2190-C  MEASURED,  THEN  ADD  THE  CURRENT  REGION  TO  THE  TOTAL  AREA  FOR  C 

2200-C  THAT  REGION.  FOUND  IS  A  FLAG  INDICATING  WHETHER  OR  NOT  AN  C 

2210-C  EXISTING  REGION  WAS  SUPPLEMENTED  OR  A  NEW  ONE  WAS  ENCOUN-  C 

2  2  2 0-C  TERED.  FOUND-1  MEANS  DATA  HAS  BEEN  STORED  IN  THE  ARRAY  AND  C 

2230-C  ANOTHER  HISTOGRAM  MUST  BE  READ.  FOUND-O  INDICATES  DATA  HAS  C 

2240-C  NOT  BEEN  STORED  AND  THE  NEXT  REGION  MUST  BE  COMPARED.  C 

2250-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


2260-C 

2270- 

STORE(I,3)-STORE(I,3)+REAL(AREA)*CAL 

2280- 

STORE (I ,4)-ST0RE(I, 4)+l .0 

2290- 

FOUND-1 

2300- 

END  IF 

2310- 

IF  (FOUND. NE.l)  THEN 

2320- 

GO  TO  100 

2330- 

END  IF 

2340- 

END  IF 

2350-C 

2360- 

FOUND-O 

2370- 

AREA-0 

2380- 

CO  TO  1 

2390-C 

2400-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2 4 1 0-C  WHEN  THE  NEXT  CARD  IN  THE  SEQUENCE  SPECIFIES  A  NEW  DAY,  THE  C 
242 0-C  PREVIOUS  DAY’S  INFORMATION  MUST  BE  SUMMARIZED.  C 
24  30-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2440*0 
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2450-  ELSE 

2460-C 

2470-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2480*0  FOR  EACH  REGION  OBSERVED  IN  A  DAY,  DETERMINE  THE  AVERAGE  C 
2490*0  AREA  (-TOTAL  AREA/NUMBER  OF  OBSERVATIONS).  C 

2500-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2510*0 

2520-  DO  200  I-l.NOREG 

2530-  STORE(I, 5) -STORE (I , 3) /STORE (1 ,4) 

2540-  200  CONTINUE 

2550-C 

2560-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
2570-C  SEARCH  THROUGH  THE  ALFA  MATRIX  UNTIL  THE  DATE  OF  THE  DAY  *  S  C 
2580-C  DATA  BEING  SUMMARIZED  MATCHES  THE  DATE  OF  THE  APPROPRIATE  C 
2590-C  P,  BO,  AND  LO.  USE  THESE  VALUES  FOR  SUMMARY  CALCULATIONS.  C 
2600-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


2610-C 

2620- 

DO  205  N-l , EPHDAY 

2630- 

DAY-INT( ALFA(N , 1 ) ) 

2640- 

IF  (DAY.EQ. YESTDY)  THEN 

2650- 

P-ALFA(N, 2)*(RAD/360) 

2660- 

B0-ALFA(N, 3  )  *  ( RAD/ 360) 

2670- 

END  IF 

2680-  205 

2690-C 

CONTINUE 

2700-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


2710-C  FOR  EACH  REGION  ENCOUNTERED,  ACCOUNT  FOR  FORESHORTENING  C 
2720-C  USING  THE  SECANT  RULE.  CALCULATE  THE  CORRECTED  AREA  (COL  7,  0 
2730-C  STORE (1,6))  AS  IF  THE  SUN  WERE  BEING  VIEWED  DIRECTLY  ABOVE  C 
2740-C  THAT  POINT.  THETA  IS  THE  ANGLE  BETWEEN  THE  CENTER  POINT  OF  C 
2750-C  THE  SUN  AND  THE  CENTER  POINT  OF  THE  OBSERVED  RECION.  C 
2760-C  IT  GIVES.  AN  INDICATION  OF  HOW  FAR  FROM  THE  CENTERLINE  OF  0 
2770-C  THE  SUN  THE  VIEWED  REGION  IS.  THE  SECANT  CORRECTION  IS  C 
2780-C  INACCURATE  FOR  LONCITUDES  GREATER  THAN  65  DECREES.  C 
2790-C  DATA  FOR  LONGITUDES  GREATER  THAN  90  DECREES  MUST  BE  C 
2800-0  DISCARDED  *  C 


28I0-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

2820-C 

2830-  DO  210  l-l , N 0 R E C 

2840-  A-ABS(STORE(I , 2)) 

2850-  B-(RAD/4)-ST0RE(I, I) 

2860-  C- ( RAD/ 4 )-B0 

2870-  THETA-ACO  S (C0S(B)*C0S(C)  +  SIN(B)*SIN(C) *COS ( A ) ) 

2880-  $TORE(I ,6)-ST0RE(I, 5  )  *  (  l .0/C0S(TUETA)) 

2890-C 

2900-CCCCCCCCCCCCCCCCCCCCCCC‘CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


2910-C  ACCOUNT  FOR  SOLAR  INCLINATION  TO  THE  ECLIPTIC  AND  ITS  VAR-  C 
292 0-C  IATION  DUE  TO  THE  EARTH’S  ROTATION  ABOUT  THE  SUN.  THIS  C 
2930-C  COORDINATE  TRANSFORMATION  GIVES  THE  LATITUDE  AND  LONGITUDE  C 
2940-C  OF  ACTIVE  RECIONS  AS  VIEWED  FROM  EARTH  IN  THE  REFERENCE  C 
2950-C  FRAME  OF  THE  ALCONQUIN  DATA  (LAT-COL  8,  LONG-COL  9).  C 


2960-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

2970-C 
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2980-  L0S-C0S(B0)*C0S(ST0RE(I, i ) ) *COS ( S TORE ( I , 2)) 

2990-  +  -SIN(B0)'SIN<ST0RE(I,1)) 

3000-  STORE (I , 9)-SIN(P)*SIN(B0)*C0S(ST0RE(I, l ) ) *C0S ( ST0RE( I , 2 ) ) 

3010-  +  +COS(P)*COS(STORE(I , 1 ) ) *S IK ( STORE ( I , 2) ) 

3020-  +  +SIN(P)*COS(BO)*SIN(STORE(I, 1)) 

3030-  Z-SIN(BO)*COS(P)*COS(STORE(l , 1 ) ) *COS ( S TORE ( I ,2)) 

3040-  +  ~SIN(P)*C0S(ST0RE(I, 1))*SIN(ST0RE(I, 2)) 

3050-  +  +C0S (P)*COS(BO)*SIN(STORE(I , l)) 

3060-  STORE(I,7)-Pl/2-ACOS(Z/SQRT(LOS**2+STORE(I,9)**2+Z**2)) 

3070*  STORE(I, 8 )-ATAN( STORE ( I ,9) /LOS) 

3080-  210  CONTINUE 

3090-C 

3100-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
3U0-C  PRINT  OUT  DAILY  SUMMARY  TABLE •  LATITUDE  AND  LONGITUDE  IN  C 
3120-C  DEGREES.  AREA  IN  MILLIONTHS  OF  THE  SOLAR  DISK.  C 

3130-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
3140-C 

3150-  PRINT*,'  ' 

3160-  PR  I  NT  *  (  **  SUMMARY  TABLE  FOR  DA  TE  :  "  ,  1 7  ,  -  IS  :  "  )  '  ,  YESTDY 

3170-  WR ITE ( * , 2  5C ) 

3180-  250  FORM AT (2X,  'REGION* ,  2X  , 'LAT* ,  4X , *  LONC '  ,  4X , 1  TOTAL  AREA*, 

3190-  +  2X, 'OBSV' ,2X, ' AVG  AREA 3X COR  AREA' ,3X, *NLAT* , 

3200-  +  4X, ,NL0NG,,3X, 'SOLRAD') 

3210-C 

3220-  DO  260  I-l.NOREG 

3230-  STORE(I, 1)-ST0RE(I, 1)*(360/RAD) 

3240-  ST0RE(I,2)*ST0RE(I, 2)*(360/RAD) 

3250-  STORECI, 7)-ST0RE(I , 7 )*(360/RAD) 

3260-  STORE ( 1 , 8 ) -STORE ( 1 , 8 ) * ( 3  60/RAD ) 

3270-  260  CONTINUE 

3280-C 

3290-  DO  230  I-l.NOREG 

3300-  WR I TE ( * , •<■ •  • •tlX,A4,2XiF7.2,lXfF7.2,lX,Fll.2,2X,F3.0# 

3310-  +  2X,F10.2,IX,F10.2,1X,F7.2,1X,F7.2,3X,F4.2) ') 

3320-  +  (RECONS (I )) , ( STORE ( I , J ) , J-l , 9 ) 

3330-  230  CONTINUE 

3340-C 

3350-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
3360-C  INITIALIZE  STORAGE  ARRAY  AND  VARIABLES  IN  PREPARATION  FOR  C 
3370-C  ANOTHER  DAY’S  DATA.  C 

3380-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCdCCCCCCCCCCCCCCCCCCCCC 
3390-C 

3400-  DO  270  I-l.NOREG 

3410-  DO  280  J-I ,9 

3420-  STORE ( I , J ) -0 . 0 

3430-  280  CONTINUE 

3440-  270  CONTINUE 

3450-C 

3460-  DO  275  1-1,14 

3470-  REGONS ( I ) -  *  • 

3480-  275  CONTINUE 

3490-C 

3500-  NOREC-O 

3510-  P-0.0 

3520-  BO-O.O 

3530-  A-0.0 

3540-  B-0.0 

3550-  C-0.0 

3560-  Z-0.0 

3570-  LOS-O.O 

3580-  TKETA-0.0 

3590-  YESTDY-YRMODA 

3600-C 
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3610-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
3620-C  SINCE  THE  FIRST  CARD  OF  A  NEW  DAY  CAUSED  THE  SUMMARY  REPORT  C 
3630-C  TO  BE  CALCULATED,  RETURN  AND  START  WITH  THE  DATA  FOR  THIS  C 
3640-C  HISTOGRAM  ON  A  NEW  DAY.  C 

3650-CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
3660-C 

3670-  GO  TO  2 

3680-  END  IF 

3690-  300  STOP 

3700-  END 

371 0-*E0R 
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APPENDIX  E 


Coordinate  Transf orma tion 

This  coordinate  transformation  algorithm  was  used  in  the 
SOON  Program  to  transpose  the  latitude  and  longitude  from  the 
SOON  coordinate  frame  to  the  Algonquin  frame  to  aid  in  direct 
comparison  of  the  data.  Xf  Y,  Z  is  the  heliocentric  frame. 
The  east-west  frame  is  x,  yf  z. 


Figure  11.  Frame  Orientations 


Rotations 

The  first  transformation  is  a  positive  or  negative  rota¬ 
tion  about  Y  by  an  angle  Bq. 
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Figure  12.  Rotation  about  Y 


The  second  transformation  is  a  positive  or  negative  rotation 


about  x  by  an  angle  P  (Ref  25:8-13). 


y  -  0 

z  0 


cos  P  sin  P  y  *  B  y 

-sin  P  cos  P  z^_  _z'_ 


therefore , 


X 

Br  y 


(27) 


where 


cos 

Bo 

0 

-sin 

Bo 

sin 

p 

sin 

B 

0 

cos  P 

sin 

P  cos  B 

0 

sin 

B 

0 

cos 

p 

-sin  P 

cos 

P  cos  B 

0, 

To  obtain  Xt  Yf  and  Z  a  spherical  to  rectangular 
coordinate  transformation  must  be  done  (Ref  34:36) 


Figure  14.  Spherical  to  Rectangular  Coordinate 

Transformation 
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S'Uch  that 


X  ■  R#  sin  ♦  cos  8 

*  R  cob  (L  )  cos  (A  ) 

•  r  r 

where 

R  *  radius  of  the  sun, 

a 

L  •  heliocentric  latitude  of  the  region  on  the 
r 

A  »  central  meridian  angle  (heliocentric)  . 
r 

Similarly, 

Y  *  R#  cos  (l»r)  sin  ( A r ) 

and 

Z  *  R#  sin  (Lr) 

Therefore,  x 

x  *  cos  R@  cos  ( Lr )  cos  (Ay)  -  sin  B0  R@  sin  (l<r)  > 

y  «  sin  P  sin  Bq  R#  cos  (Lr)  coa  (Ar) 

+  cos  P  R^  cos  (L_)  sin  (A  ) 
e  r  * 

♦sin  P  cos  B^  R_  sin  (L  )  , 

o  •  r 

and 

z  ■  sin  B  cos  P  R^  cos  (L  )  cos  (A  ) 

o  •  r  r 


(29) 


sun,  and 


(30) 


(31) 


(32) 


(33) 


The  data  was  unitized  to  solar  radii  by  setting  R  ®1 .0  . 

© 

To  go  back  to  latitude  and  central  meridian  angle  in  the 
east-vest  frame,  a  rectangular  to  spherical  transformation  is 
performed : 

tan  >r  -  y  /  x  (35) 

then  *  , 

Xr  -  tan_1(y  /  x) ,  x  7  0  (36) 

and 

cos  (90  -  i  )  -  *  /  /x7+yZ+z7  (37) 

therefore 

lr  -  90  -  cos_1[z  /  /x2+y7+z7  ]  (38) 

where  \ 

Xr  *  east-west  frame  central  meridian  angle  (longitude) 

and 


1  *  east-west  frame  latitude* 

r 

In  the  SOON  Program 

x  -  LOS  (line  of  site),  and 
y  -  STORE  (1,9). 
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APPENDIX  P 


Da ta  Analysis  Programs 

The  follvoing  programs,  ANAL1  and  ANAL3f  were  used  to 
analyse  the  data  for  this  project.  ANAL1  provides  a  scatter- 
gram  of  the  data  while  ANAL3  gives  a  more  complete  breakdown 
of  the  results  of  the  analysis  including  the  residuals  used  to 
determine  the  adequacy  of  the  model.  These  basic  programs, 
with  minor  modifications,  were  used  to  test  the  data  after 
removing  some  of  the  data  points.  ANAL3  was  used  to  find  more 
significant  digits  for  the  predictive  equations  by  multiplying 
each  NEWSCN  and  NEWAREA  value  by  1000  on  one  run. 
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1 00-CAP, CM60000,T100, 10120.  T820740 , PUZ , *172 
1 10- ATTACH , SPSS, XD-A FIT. 

120-SPSS. 

130-*E0R 

140-RUN  NAME  CALSCN  VS  S00NAREA  N0V82  ANAL1 

150-VARIABLE  LIST  OATE, CALSCN, SOONAREA 

160-VAR  LABELS  CALSCN  CALIBRATED  OTTAWA  SCAN  VALUES/ 

170-  SOONAREA  CORRECTED  AREA  FROM  SOON  HISTOGRAMS 

180-INPUT  FORMAT  FREEFIELD 

190-N  OF  CASES  142 

200-COMPUTE  NEWSCN-CALSCN/ SOONAREA 

2 10-C0MPUTE  NEWAREA-1 /SOONAREA 

220-PRINT  FORMAT  CALSCN  (2), SOONAREA  (2),NEWSCN  (6) .NEWAREA  (6) 
230-LIST  CASES  CASES-142/VARIABLES-DATE, CALSCN, SOONAREA, 

240-  NEWSCN, NEWAREA 

2 50-SCATTERGRAM  NEWSCN  WITH  NEWAREA 

260-STATISTICS  ALL 

270-READ  INPUT  DATA 

280-780501  4.82  13686.81 

290-780504  4.70  26714.89 

300- ... 

310-FINISH 
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*100-CAP,CM60000,T100, 10120.  T820740, PUZ , 4172 
UO-ATTACH,  SPSS, ID-AFIT. 

120-SPSS. 

130-*E0R 

140-RUN  NAME  CALSCN  VS  S00NAREA  N0V82  ANAL3 

ISO-VARIABLE' LIST  DATE , CALSCN, SOONAREA 

160-VAR  LABELS  CALSCN  CALIBRATED  OTTAWA  SCAN  VALUES/ 

170-  SOONAREA  CORRECTED  AREA  FROM  SOON  HISTOGRAMS 

180-INPUT  FORMAT  FREEFIELD 

190-N  OF  CASES  142 

200-C0MPUTE  NEWSCN-CALSCN/ SOONAREA 

210-COMPUTE  NEWAREA-I /SOONAREA 

220-PRINT  FORMAT  CALSCN  (2), SOONAREA  (2),NEWSCN  ( 6 ) , NEUAREA  (6) 
230-LIST  CASES  CASES-1 42/VARIABLES-DATE , CALSCN, SOONAREA, 

240-  NEW SCN , NEW AREA 

2S0-REGRESSI0N  METHOD-FORWARD/ 

260-  VARIABLES-NEWSCN, NEW AREA/ 

270-  REGRESSION -NEW SCN  WITH  NEWAREA/ 

280-  RESIDUALS/ 

290-STATISTICS  ALL 

300-READ  INPUT  DATA 
310-780301  4.82  13686.81 
320-780304  4.70  26714.89 
330-.. . 

340-FINISH 
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